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SIZES AND TEMPERATURE PROFILES OF QUASAR ACCRETION DISKS 
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ABSTRACT 

Microlensing perturbations to the flux ratios of gravitationally lensed quasar images can vary with wavelength 
because of the chromatic dependence of the accretion disk's apparent size. Multiwavelength observations of 
microlensed quasars can thus constrain the temperature profiles of their accretion disks, a fundamental test of 
an important astrophysical process which is not currently possible using any other method. We present single- 
epoch broadband flux ratios for 12 quadruply lensed quasars in 8 bands ranging from 0.36 to 2.2 /im, as well as 
Chandra 0.5 - 8 ke V flux ratios for five of them. We combine the optical/IR and X-ray ratios, together with X- 
ray ratios from the literature, using a Bayesian approach to constrain the half-light radii of the quasars in each 
filter. Comparing the overall disk sizes and wavelength slopes to those predicted by the standard thin accretion 
disk model, we find that on average the disks are larger than predicted by nearly an order of magnitude, with 
sizes that grow with wavelength with an average slope of ^0.2 rather than the slope of 4/3 predicted by the 
standard thin disk theory. Though the error bars on the slope are large for individual quasars, the large sample 
size lends weight to the overall result. Our results present severe difficulties for a standard thin accretion disk 
as the main source of UV/optical radiation from quasars. 
Subject headings: accretion, accretion disks - gravitational lensing: strong - quasars: general 



1. INTRODUCTION 

Microlensing of the images of gravitationally lensed 
quasars by stars in the foreground galaxies causes time- 
variable anomalies in their flux ratios (see review by Wamb- 
sganss 2006). In-depth study of these anomalies offers a 
unique glimpse into the properties of both the deflector and 
the quasar itself. Because the amplitude of microlensing 
anomalies depends on the angular size of the source, it can 
be used to estimate the size of quasar emission regions at the 
level of microarcseconds, far beyond current resolution limits 
(e.g., Rauch & Blandford 1991 ). The temperature structure of 
quasar accretion disks on these same scales causes chromatic 
microlensing, where blue light from the inner regions is more 
strongly microlensed than red light from farther out ( Wamb- 
|sganss & Paczy nski 1991]). Thus, multiwavelength observa- 
tions of quasar microlensing have a unique capability to mea- 
sure the structure of accretion disks and put direct constraints 
on accretion models. 

Several recent studies have made use of microlensing to ex- 
plore the size sc ale of the accretion d isk, which emits most of 
the op tical light. Pooley et al.l (2007) fo llow up on earlier re- 
sults (Blac kburne et al.|20 06; Pool ey et al|2006| ) to show that 
the flux ratio anomalies caused by microlensing are in general 
stronger in X-rays than at optical wavelengths, and infer that 
the source of optical light is larger than the X-ray source, and 
larger on av erage than exp ected fo r a standard thin accretion 
disk model (Shakura & Sunyaev 1973 INovikov & Thorne 



1973 1. Morgan et al. (2"0l0) estimate disk sizes for several 
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lensed quasars using time-series observations of microlens- 
ing variability, rather than single-epoch flux ratio anomalies. 
They find only marginally better agreement with thin disk 
models, and note that flux-based size predictions are system- 
atically too small. 

Several other studies move beyond single-wavelength size 
estimates to estimate the dependence of the source size on 
wavel ength and thus constrain the temper ature p rofile of the 
disk. |Bate et al I ( [2008] ) and |Floyd et al.] d2009| place upper 
limits on the radii oFaccretion disks using single-epoch broad- 
band flux ratio measurements; they also put weak constraints 
on the temperature profile slope. Mosquera e t al.| ( |2010| ) 
make size and slope estimates using narrowband photome- 
try. Other studies combine multiwa velength observations with 
time-series light curves ( Anguita et al. 2008 ; Poindexter et al. 
2008 ; Mosquera et al. 2009; Poindexter & Kochanek 2010a), 
or even employ spectroscopic monitoring (Eigenbrod et al. 
2008). Additionally, efforts have been made (Morgan et al. 
2008, Dai et al. 2010) to estimate the size of the source of 
the non-thermal X-ray flux by combining Chandra measure- 
ments with optical light curves. Because of the effort involved 
in these observations, each of these studies focuses on an in- 
dividual lensed quasar. 

The effects of microlensing also depend on properties of 
the foreground galaxy — the mass function of the stars, and 
the surface density of clumpy stellar matter as a fraction of 
the total (smooth and clumpy) surface density. The amplitude 
of microlensing perturbations decreases for very small stel- 
lar fractions, and also decreases, perhaps counterintuitively, 
for large stellar fractions (Schechter & Wambsganss 2002). 
This has led to results that effectively rule out the 100% stars 
case as well as the 100% smooth matt er case, with a best- 
fit stellar mass fraction of ^10% (Schechter & Wambsganss 
2004; Pooley et al. 2009). The mass function of the microlens 
stars changes the variab ility characte r istics o f microlensing 
light curves; Poindexter & Kochanek (2010b) exploit this to 
estimate the mean stellar mass in the lensing galaxy of the 
quadruple lens Q 2237 + 0305. 
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With the goal of measuring accretion disk sizes and tem- 
perature profiles, we present single-epoch observations of the 
flux ratios of a sample of 12 lensed quasars in eight broadband 
optical and infrared (IR) filters in the "big blue bump" region 
of the spectrum. The low cost of our single-epoch observing 
strategy relative to long monitoring campaigns allows us to 
cover a factor of six in wavelength, while roughly tripling the 
number of lensed quasars with observations of this kind. We 
also present X-ray flux ratios from the Chandra X-ray Obser- 
vatory for 5 of the quasars. Comparing these flux ratios, as 
well as X-ray ratios from the literature for six of the seven re- 
maining lenses, to those predicted by smooth lens models, we 
use a Bayesian analysis method to determine posterior proba- 
bility distributions for the sizes of the quasar emission regions 
in each filter. This allows us to estimate the overall size scales 
of the accretion disks, as well as the scaling of the size with 
black hole mass and wavelength. As we will see, the rela- 
tively large size of our sample is crucial for drawing conclu- 
sions about accretion disk structure. 

Our sample is entirely made up of quadruply lensed 
quasars; they have the advantage of having three flux ratios 
apiece (rather than one for double lenses) and are easier to 

model. They are as follows: 

HE 0230 -21 30 dWisotzM et al.|1999) , 
MG J0414+053 4 (|Hewitt et al.|1992) , 
HE 0435- 1223 ((Wisotzki et al. 2002), 
RX J09 1 1 + 055 1 ( Bade e tal.|1997), 
SDSS J0924 + 02 19 dlnada et al.|2003), 
HE 1113-0641 (Blackbu rne~etal.|2008| ), 
PG 1 1 1 5 + 080 ( Weyman neTal.|1980| l, 
RXJ1131-1231 (Sluseetal. 2003), 
SDSSJ1138 + 0314(Eigenbrodetal. 2006), 
SDSS J1330+ 1810 (Oguri et al. 2008), 
WFIJ2026-4536 (Morgan et al. 2004), and 
WFIJ2033-4723 (Morg an et al.pOOfy 
For brevity, we will usually exclude the declination from their 
names from now on. Following Pooley et al. (2007), we cat- 
egorize the quasar images by their macro-magnification and 
parity: highly magnified minima and saddle points of the light 
travel time surface are designated HM and HS, while their 
less-magnified counterparts are LM and LS, respectively. 

The variation of accretion disk size with black hole mass 
and wavelength provides direct constraints on disk models. 
The standard thin disk model predicts a profile for the effec- 
tive temperature that falls approximately as the j3 = 3/4 power 
of the radius: 
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where Mbh is the black hole mass, m p is the proton mass, 
o~b is the Stefan-Boltzmann constant, and a T is the Thomson 
scattering cross-section. The dimensionless factors / E( y and 
■q are the ratios of the quasar's bolometric luminosity to its 
Eddington luminosity and to the accretion rate Mc 2 , respec- 
tively. The function g is defined as g(x) = 1 -jc 1 / 2 , and r m is 
the innermost radius of the accretion disk. As we will see, at 
optical wavelengths the disk size is much larger than r m , so 
for simplicity we set g to unity. Likewise, the general rela- 
tivistic corrections of Noviko v"& Thorne] ( |1973| l are small at 
the radii of interest. By integrating the specific flux at fixed 
wavelength over the entire multi-temperature disk, we calcu- 
late the half-light radius — that is, the radius at which half of 



the light at a given wavelength is emitted — predicted by the 
thin disk temperature profile: 
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In this equation, h P is Planck's constant, and i is the unknown 
inclination angle. We include the factor of (cos i) 1 ' 2 to account 
for the inclination of the disk to our line of sight, and set it to a 
reasonable average value of (1/2) 1 / 2 in the second line. For a 
thin disk, the half-light radius scales with the rest wavelength 
as A 4 / 3 ; for a more general disk temperature profile 7^g oc r~@, 
the half-light radius will depend on A 1///3 . 

In Sections 13] and [3] we describe the X-ray and optical/IR 
observations and flux ratio measurements. In Section [4] we 
estimate the systematic uncertainties in the measured flux ra- 
tios. We describe our lens models in Section [5] and in Sec- 
tion|6]we describe our Bayesian analysis method. We discuss 
our results and give concluding remarks in Sections [7] and M 
Throughout this work, we calculate distances and time delays 
using a geometrically flat universe with Q,m = 0.3, J7a = 0.7, 
and H = 70 km s" 1 Mpc" 1 . 

2. X-RAY FLUX RATIOS 

The X-ray data in this survey are from the Advanced CCD 
Imaging Spectrometer (ACIS) instrument aboard the Chandra 
X-ray Observatory. For five of the lenses among our sample 
(specifically HE0230, MGJ0414, RXJ0911, SDSSJ0924, 
and WFIJ2033) we adopt the X-ray flux ratios reported by 
Pooley et aT] ( |2007) . For PG 1 1 15 we adopt the X-ray flux 
ratio s of ""th e 2008 Janua ry 31 observation (ObsID 7757) 
from |Pooley et al. (2009), as they are more contemporane- 
ous with our optical/IR observations. For the same reason, 
for RX Jl 131 we use observations from the Chandra archive. 
Specifically, we use the 2007 February 13 observation (ObsID 
7787) for comparison with optical data and the 2007 April 16 
observation (ObsID 7789) for comparison with IR data (see 
Section |3J. The high level of X-r ay microlensing variability 
in these two systems in particular (Chart aset al.|20 09 , Pooley 
et al.|20 09) underscores the importance of making our multi- 
wavele ngth observations as contemporaneous as possible (see 
Section [4~2| for a detailed discussion). 

The remaining lenses in our sample, excluding 
SDSSJ1330, which has not been observed in X-rays, 
were observed using ACIS between 2006 December 17 and 
2007 June 2 8. The data were re duced using the procedure 
described by |Pooley et al] ( |2009[ ). We produced 0.5-8 keV 
images of each lens with a resolution of 0."0246 pixel" 1 , and 
fit four Gaussian components and a constant background 
to each image. The background level was fixed at a level 
determined from a nearby source-free region, and the relative 
positions of the four Gaussians were fixed to the Hubble 
Space Telescop e (HST) positions provided by the CASTLES 
database)^] (Falco et al. 2001). The width of the Gaussian 
point-spread function (PSF) was allowed to vary for each 
lens, but was constrained to be the same for all four quasar 
components. The fits were performed using |Cash| ( |1979| l 



5 See http : //www. c fa. harvard. edu /castles 
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Quasar 



ObsID c Date 


Exposure (s) 


HE 0435 -1223 


7761 2006 Dec 17 


10130 


HE 1113-0641 




7760 2007 Jan 28 


15 180 


RXJ1131-1231 d 




7787 2007 Feb 13 


5 190 


7789 2007 Apr 16 


5 190 


SDSSJ1138 + 0314 




7759 2007 Feb 13 


19 080 


WFIJ2026-4536 





Table 1 

X-ray Fluxes and Flux Ratios 

Image Flux Ratios" 



7758 2007 Jun 28 



HS/LM 
B/A 



D/A 

0.78!! 



HM/LM 
C/A 



B/A 



LS/LM 
D/A 



C/A 

0.20^: 



10 170 



A/C 



D/C 

1 3 +0 - 6 

A2/B 
2.0!!' 



B/C 



A/C 

3 2 +L0 
J.A-i.o 

Al/B 
5.81}? 



D/C 



B/C 

1 o +0 - 4 

C/B 
0.4O2'?5 



LMUnabs. F .5-8keV 
(10~ 14 ergcnT 2 s~') 

Fa 



Fc 

45.34!"' 



1 6 +09 
5.3!!? 



Note. — We only report new measurements here; see references in the text for other flux ratios. 
11 HS: highly magnified saddle point; HM: highly magnified minimum; LS: less magnified saddle point; LM: less 
magnified minimum. 

b The unabsorbed flux of the LM image is computed from the best fit power-law model described in Sectional 
c The observation identifier of the Chandra dataset. 
d The data for RXJ1 131- 1231 come from the Chandra archive. 



statistics, which are appropriate when small numbers of 
photons have been collected. In the cases of HE 1113, 
SDSSJ1138, and WFIJ2026, the least-squares optimization 
settled on a wider PSF than for any of the other observations 
in our previous work. Given the low numbers of total counts 
in these observations and the small separations of the images, 
we fixed the width of the PSF in these cases to the average 
of the other X-ray PSFs. The flux in image A2 of WFIJ2026 
was very sensitive to the PSF, so we allocated a generous 1 .5 
mag of uncertainty to its flux ratio, added in q uadrature to the 
other sources of uncertainty (see Section [42] ). The other flux 
ratios were not sensitive to the PSF size. Table [TJ summarizes 
the Chandra observations and the measured flux ratios. 

3. OPTICAL AND INFRARED PHOTOMETRY 

Between 2007 February and 2008 May, we undertook an 
optical observing campaign to obtain multi-band, contempo- 
raneous images of our sample of lenses. At IR wavelengths, 
we used the J, H, and K s filters with Persson's Auxiliary Nas- 
myth Infrared Camera (PANIC), on the 6.5 m Baade telescope 
at Las Campanas Observatory. In the optical, we used the 
Sloan u'g'r'i'z' filters with the Raymon and Beverly Sackler 
Magellan Instant Camera (MagIC), which was on the neigh- 
boring Clay telescope when we began our campaign, but was 
moved to Baade while the observations were still underway. 
This change of location allowed us to use both instruments 
during the same observing run, and even during the same 
night. MagIC and PANIC have fields of view 2f4 and 2.'1 on 
a side, respectively, large enough to include stars appropriate 
for use as PSF templates. The instruments' pixel scales are 
0"069 and 0."125 pixels, which more than adequately sample 
the PSF. The details of the observations are listed in Table [2] 

During our 2007 September observing run, MagIC was of- 
fline for an upgrade, so we instead used the Inamori Magellan 
Areal Camera System (IMACS) in its imaging mode. With its 
f/4 camera, the instrument has a pixel scale of 0" 1 1 and a 
15. '5 field of view, but in order to reduce the readout time we 
used only a subraster 2.'2 on a side. We used this instrument 
for our observations of HE 0230, MG J0414, and HE 0435. 



The images were corrected for bias and dark current, 
flattened, and cleaned of cosmic rays using standard tech- 
niques. Where multiple exposures were obtained, they were 
stacked, resulting in a single image per filter per lensed quasar. 
Postage-stamp /'-band images of the 12 lenses in our sample 
are provided in Figure [T] together with the residuals from the 
PSF fitting described irTSection[3] 

Because of the small separations between the quasar images 
in our sample, we only observed under the best atmospheric 
seeing conditions, typically below 0"7. The excellent image 
quality at the Magellan telescopes was crucial to the success 
of our survey. The i! - and /-band seeing is reported in Table[2j 
which outlines the optical/IR observations obtained for this 
work. The seeing in the other filters was roughly consistent 
with the rule of thumb: FWHM ex A" 1/5 . 

Because we are primarily interested in the flux ratios of 
lensed quasars (indeed, we cannot measure absolute magni- 
fications), we did not, in general, obtain images of standard 
stars for the purpose of calibrating our photometry. In the 
case of SDSS J1330, however, we did observe standard stars 
in the optical bands, so in Section|3]we report calibrated pho- 
tometry for this lens. We also estimate the photometric zero 
point in the J and K s bands, based on images taken the same 
night which contain several calibra ted source s from the Two- 
Micron All Sky Survey (2MASS; |Skrutskie et al.|2006) >. 

Despite the high quality of our data, careful PSF subtrac- 
tion was necessary to disentangle the compact clumps of point 
sources and lens galaxies that make up our sample. For each 
image, we used a nonlinear least-squares fitting algorithm to 
simultaneously fit the positions and relative fluxes of the four 
quasar components and the lens galaxy. 

For lenses with well-separated components, we performed 
a simultaneous fit to the images in all eight bands, with the rel- 
ative positions of the quasar components and the lens galaxy 
allowed to vary, but constrained to be consistent across all 
eight filters. The flux ratios of the quasar components were 
free to vary in each filter independently. The best-fit positions 
were consistent in every case with the HST positions from the 
CASTLES survey. In cases where merging pairs or triplets 
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Table 2 










Optical 


Observations 






Quasar 


Date 


Instrument 


Filters 


Exposures (s) 


Seeing 


HE 0230 


2007 Sep 16 


JMACS 


u'g'r'i'z' 


720; 240; 240; 240; 480 


0"43 




2007 Jul 29 


PANIC 


JHK, 


540; 540; 540 


0"64 


MGJ0414 


2007 Sep 21 


IMACS 


r'i'z' 


360; 240; 480 


0"67 




2007 Sep 22 


PANIC 


JHK, 


810;720;720 


0"76 


HE 0435 


2007 Sep 16 


IMACS 


u'g'r'i'z' 


720; 240; 240; 270; 480 


0"67 




2007 Sep 22 


PANIC 


JHK, 


810; 1080; 1080 


1"06 


RXJ0911 


2007 Feb 13 


MagIC 


g'r'i'z' 


240; 240; 300; 480 


0"70 




2007 Feb 14 


MagIC 


u' 


720 


0"63 




2007 Apr 7 


PANIC 


JHK, 


540;405;405 


O'.'H 


SDSS J0924 


2007 Feb 13 


MagIC 


u' g'r'i'z' 


720; 240; 240; 180; 480 


0."61 




2007 Apr 7 


PANIC 


JHK, 


540; 540; 540 


0"64 


HE 11 13 


2007 Feb 14 


MagIC 


g'r'i'z' 


240;180;120;240 


0"49 




2007 Apr 7 


PANIC 


JHK, 


450; 450; 450 


0"45 


PG1115 


2008 Feb 1 


MagIC 


g'r'i'z' 


120; 120; 150; 240 


0"55 




2008 Feb 2 


MagIC 


u' 


720 


0."87 




2008 Feb 4 


PANIC 


JHK, 


540; 540; 540 


0"56 


RXJ1131 


2007 Feb 13 


MagIC 


u'g'r'i'z' 


480; 120; 120; 120; 240 


0."65 




2007 Apr 6 


PANIC 


JHK, 


486;405;405 


0"76 


SDSSJ1138 


2007 Feb 13 


MagIC 


u'g'r'i'z' 


1080;480;480;360;960 


0"75 




2007 Apr 8 


PANIC 


JHK, 


540; 540; 540 


0"47 


SDSSJ1330 


2008 Feb 1 


MagIC 


u'g'r'i'z' 


360;120;120;150;240 


0"53 




2008 Feb 3 


PANIC 


JHK, 


540; 540; 540 


0"68 


WFIJ2026 


2008 May 13 


MagIC 


u'g'r'i'z' 


720; 240; 240; 360; 480 


0"54 




2008 May 12 


PANIC 


JHK, 


135;135;90 


0"43 


WFIJ2033 


2007 Jun 15 


MagIC 


u'g'r'i'z' 


720; 240; 120; 300; 480 


0"88 




2007 Jul 4 


PANIC 


JHK, 


135;135;90 


Of'51 



of images would have caused strong correlations between po- 
sitions and flux ratios, we fixed the relative positions of the 
quasar components and the lensing galaxy to the CASTLES 
positions. We modeled the quasar images using nearby rela- 
tively bright stars as empirical PSF templates. 

In the r' and redder bands, we did not use least-squares min- 
imization to determine the flux ratios of RX J 1131. Because 
of the bright Einstein ring in this system, this strategy would 
have over-subtracted the quasar components. Instead we fixed 
all of the quasar fluxes to values that resulted in residuals that 
looked like an unbroken Einstein ring. Although it was com- 
pletely ad hoc, we think that this "chi by eye" techniq ue g ave 
flux ratios less affected by systematic errors. Section |4~T| de- 



scribes our method for estimating the uncertainty of these flux 
ratios. 

In all cases we modeled the lensing galaxy as a two- 
dimensional pseudo-Gaussian function (Schechter et al. 
1993), with a full-width at half-max as a fixed parameter. The 



widths were chosen using a trial-and-error technique, exam- 
ining the residuals by eye. (Because of the frequent presence 
of a faint Einstein ring due to the quasar's host galaxy, least- 
squares minimization often overestimated the width.) Be- 
cause the lens galaxy was nearly always faint compared to 
the quasar images, the goodness of the fit was insensitive to 
this approximation. The galaxies were mostly round, but for 
the few exceptions we fixed the axis ratio and position angle 
by hand, as we did the width. 

The fitted fluxes of the lens galaxies are uncalibrated, and 
therefore of little interest; we do not report them. SDSS J 1 330 
is the exception; because it was discovered relatively recently, 
and because observing conditions were photometric on the 
night it was observed, we observed standard stars in the opti- 
cal bands, for the purpose of calibrating its fluxes. The cali- 
brated magnitudes of this lens galaxy are reported in Table HI 
The optical calibration relies on aperture photometry of the 
standard stars SA 107-351 (for u') and G163-50 (for g'r'i'z'), 



and places the fluxes on the Sloan u'g'r'i'z' photometric sys- 
tem, which closely ap proximates the monochromatic AB sys- 
tem ( |Smith et al.|2002) . In the H and K s bands, we made use 
of a field containing several 2MASS sources for photometric 
calibration; we performed aperture photometry on them and 
compared their measured magnitudes to those reported in the 
NAS A/IPAC Infrared Science Archive^] to obtain photomet- 
ric zero points. This places our calibrated / and K s fluxes on 
the 2MASS photometric system ( |Cohen et al.|2003) . In the 
optical bands, used the atmospheric extinction coefficients of 
Smith et al. ( |2007| l to correct for differences in airmass be- 
tween SDSSJ1330 and the standard stars; these corrections 
were all less than 0.04 mag. 

The lack of HST positions for SDSS J 1330, combined with 
the relative brightness and elliptical shape of the lensing 
galaxy, made its fit a special challenge. In order to get the 
shape of the galaxy right, we adopted an iterative approach. 
First, we subtracted a rough model of the lens galaxy, fixing 
its parameters to reasonable guesses. We then fit the posi- 
tions and fluxes of the quasar components, subtracting them 
from the image. We added back our rough galaxy model to 
the residual image and fit a single pseudo-Gaussian to the 
resulting galaxy image, with all parameters allowed to vary. 
The best-fit FWHM of the galaxy was 1 "0 along the major 
axis; the axis ratio and position angle were 0.67 and 24?5 
east of north. Finally, we fit the original image again, fixing 
the galaxy shape parameters to these best-fit values, in order to 
find the quasar positi ons and flu xes. The results are consistent 
with those of Oguri et al. ( 2008 1, but have higher precision; 
they are reported in Tables|3](astrometry) andHl(photometry). 

For each lens, the transformation between pixel coordinates 
and sky coordinates was determined by registering our i'- and 
/-band images to the USNO-B astrometric catalog. On aver- 
age, there were ~ 10-15 objects for each fit. The pixel scale 



http://irsa.ipac. caltech.edu 
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Figure 1. Postage-stamp images and residuals of lenses in the i band. From 
top left, they are HE0230, MGJ0414. HE0435, RXJ0911, SDSSJ0924, 
HE1113, PG1115, RXJ1131, SDSSJ1138, SDSSJ1330, WFIJ2026, and 
WFIJ2033. Each frame is 5" on a side, with north up and east left. The 
gray-scale stretch on the residual images is set to ±20<r of the sky level. For 
RXJ1 131, we show the residuals of the "chi-by-eye" fit, where the flux ratios 
are set by hand to account for the bright Einstein ring. 



and rotation angle thus determined for i! and J were assumed 
to apply to the remaining MaglC/IMACS data and PANIC 
data for that lens, respectively. We tested this assumption for 
a few cases by fitting images in other bands to the catalog and 
found it to be justified. 

Relative photometry for all the lenses in our sample except 
SDSS J1330 is listed in Table We also list the calculated 
rms deviation of the magnitudes in each filter from those pre- 
dicted by the lens models. This number is an indication of 
how anomalous the flux ratios are. For comparison with the 
optical flux ratios, we also list the X-ray ratios, in magnitudes. 

4. UNCERTAINTY ESTIMATION 

Aside from statistical measurement uncertainties in the flux 
ratios, systematic errors arise from several sources. At optical 
and IR wavelengths, quasar emission lines, quasar variability, 
or contamination from a lens galaxy or an Einstein ring may 
contribute. In X-rays, quasar variability is a factor, and delays 
between X-ray and optical/IR observations combine with mi- 
crolensing variability to contribute additional systematic un- 
certainty. 



Table 3 

Relative Astrometry for SDSS J1330+ 1810 

Component x y 



A 


-1.26±0.01 


-1.17±0.01 


B 


-0.84 ±0.01 


-1.19±0.01 


C 


= 


= 


D 


-1.49 ±0.01 


0.44 ±0.01 


G 


-1.04 ±0.02 


-0.20 ±0.02 



Note. — All positions are in arcseconds. 
Positive x and y point west and north, respec- 
tively. 



4.1. Optical Uncertainties 
Quasar emission lines are thought to come fro m a region 



& 



Schneider 



too large to be s trongly affected by microlensin^F 

Wambsganss 1990). Therefore the presence of emission 
line flux in our broadband measurements will cause errors if 
we assume we are measuring only the continuum from the 
accretion disk. The strength and effect of these errors is very 
difficult to predict. We allocated a 0.05 mag systematic un- 
certainty in our flux ratios for filters into which one of the 
following lines has been redshifted: Civ, Cm], Mgll, H/3, or 
Ha. (This is very roughly the percentage of the broadband 
flux taken up by one of these lines, in general.) Occasionally, 
Lya falls in a filter, or there are two emission lines present; in 
these cases, we allocated 0. 1 mag of uncertainty. 

The multiple images of a strongly lensed quasar arrive with 
relative delays of hours to weeks because of the different paths 
taken by their light. Quasar variability can conspire with these 
time delays to mimic flux ratio anomalies. We do not ex- 
pect this to be a very strong effect, because quad lenses tend 
to have short time delays, especially between close pairs and 
triplets where the most interesting anomalies tend to happen. 
In order to quantify the effect, we extrapolated the empirical 
quasar variability structure function in Figure 18 of |de Vries] 
et al. (2005 i using a power law: logS(r) = 0.8 + 0.651ogr. 
This gave us an estimate for the standard deviation S of the 
quasar brightness in magnitudes as a function of time delay r 
in years. 

Time delays have been measured for five of the lenses in our 
sample. For HE0435, RXJ0911, PG1115, RXJ1131, and 
WFI J2033, we used, respectively, the time delays reported by 



Kochanek et al. (2006), Hjorth et al. (2002), Barkana ( 1997), 
Morgan et al. (2006), and Vuissoz et al.|( |2008| ). For the re- 
maining lenses we used the time delays predicted by our lens 
models (see Section B). The lens redshifts of HE 1113 and 
WFIJ2026 are unknown, so we estimated them both to be 



Zl ~ 0.7 using the method of Kochanek ( 1992 1. The resulting 
quasar variability uncertainties were only significant for a few 
quasar images (e.g., image D in RXJ1131 or RX J0911). 

Finally, the formal measurement uncertainties for a few flux 
ratios were clear underestimates, or (in one case) simply did 
not exist. HE 0230 provides an example, as its component D 
is blended in our images with a bright companion to the main 
lens galaxy only 0."4 away. For this flux ratio, we allocated 
an extra error equal to the change in measured brightness if 
25% of the galaxy light were attributed to the quasar image. 
This ranged from 0. 1 mag in g' to 1 .6 mag in K s , The flux ra- 
tios of RX J1131 have no formal uncertainties because of the 

7 Some degree of broad line emission microlensing has been observed, 
however (e.g., Keeton et al. 2006 1. 
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Table 4 

Photometry for SDSS J1330+ 1810 



Filter 


A-C 


B-C 


C 


D-C 


G 


u' 


-0.75 ±0.01 ±0.05 


+0.27 ± 0.02 ± 0.05 


19.69 ±0.06 ±0.01 


+1.70 ±0.02 ±0.05 




s' 


-0.79 ±0.01 ±0.05 


+0.17 ±0.01 ±0.05 


20.99 ±0.06 ±0.01 


+1.78 ±0.02 ±0.05 


25.08±0.11 


r' 


-0.86 ±0.01 ±0.05 


-O.18±0.01±0.05 


20.57 ±0.06 ±0.01 


+1.63 ±0.02 ±0.05 


23.48±0.11 


i' 


-0.90 ±0.01 ±0.00 


-0.13±0.01±0.00 


21. 07 ±0.06 ±0.01 


+1.71±0.02±0.01 


23.26±0.10 


z' 


-0.88 ±0.01 ±0.00 


-0.12±0.02±0.00 


21.55 ±0.06± 0.01 


+1.95 ±0.05 ±0.01 


23.46 ±0.10 


J 


-0.85 ±0.02 ±0.05 


-0.18 ±0.03 ±0.05 


18.88±0.06±0.01 


+1.33 ±0.06 ±0.05 


18.94±0.10 


H 


-0.84 ±0.03 ±0.05 


-0.24 ±0.03 ±0.05 




+1.54 ±0.09 ±0.05 




K s 


-0.89 ±0.04 ±0.00 


-0.36 ±0.04 ±0.00 


17.78±0.06±0.01 


+2.29 ±0.30 ±0.01 


17.13±0.10 



Note. — Optical magnitudes for image C and lens galaxy G are calibrated to the Sloan u'g'r'i'z' 
system, while J and K s magnitudes are calibrated to the 2MASS system. Formal statistical uncertainties 
are reported first, followed by estimated systematic uncertainties. No systematics are reported for the lens 
galaxy G. 



ad hoc method we used to estimate them (see Section [3]). We 
conservatively estimated their uncertainties as half the differ- 
ence between our estimates and those of the poorly perform- 
ing chi-square fit. In addition, we allocated 0.5 mag of extra 
uncertainty to images C and D of SDSS J0924 in the u' band, 
which were only marginally detected, and 0.05 mag to images 
Al and A2 of WFI J2026 in the H band, whose flux ratios de- 
pended slightly on our choice of PSF comparison star. 

Table[5jshows the formal statistical uncertainty of each flux 
ratio, followed by the systematic uncertainty, defined as the 
quadratic sum of that due to emission lines, intrinsic variabil- 
ity, and other systematics, as described above. 



4.2. X-ray Uncertainties 

Like the optical flux ratios, the X-ray ratios had measure- 
ment uncertainties due to measurement noise and the blending 
of close pairs; both errors were generally larger because the 
X-ray observations have fewer photons and a broader PSF. 
There were also contributions from intrinsic quasar variabil- 
ity, as in the optical case. The X-ray ratios had no system- 
atic errors due to emission lines, but they did have errors due 
to microlensing variability because they were not generally 
measured contemporaneously with the other wavelengths. 

Crucial to our analysis is the assumption that the arrange- 
ment of the source and the microlenses is the same for all 
wavelengths. But when observations are not contemporane- 
ous, the source and the microlenses have the opportunity to 
reconfigure themselves. To estimate the magnitude of this ef- 
fect, we again used a structure function; this time it was not 
empirical, but was derived from microlensing magnific ation 
patterns (see our description of these patterns in Section |oTT| i. 
For each image, we chose the appropriate patt ern ( with the 
stellar mass fraction k* set to 0.1; see Section |6.1| l and ran 
1000 tracks across it, in random directions, measuring the 
variance in the logarithm of the magnification as a function 
of distance moved. 

The conversion to a structure function (with a time delay on 
the abscissa instead of a distance) required an estimate of the 
transverse speed of the source relative to the lens. We added 
four velocity components in quadrature: the velocity disper- 
sion of the stars in the lens galaxy (as estimated by our lens 
model), the tangential component of the velocity of the Sun 
relative to the rest frame of the cosmic microwave background 



(CMB), and the peculiar velocities of the qua sar and the len s 
galaxy (as estimated using Equation (14.10) of Peebles [1980[ l: 



0"L ^os 



(l+z L )D 0L 



+ v-cMBSina 



£>L S y 



+2 



/fe)V +: 



d+^) 3/2 /(0)y 



DolJ 

gpec f(Zl) Dos \ 

(1+Zl) 3 / 2 /(0)£>oJ 



(3) 



Here zs and z L are the redshifts of the source and the lens, re- 
spectively, and Dol, Aas, and Dls are angular diameter dis- 
tances from observer to lens, observer to source, and lens to 
source. We used these distances to project all velocities to the 
source plane. The angle a is measured between the Sun's ve- 
locity with respect to the CMB rest frame an d the line of sight 
to the lens, and we s et vcmb to 370 km s _1 ( |Lineweaver et al. 
|1996| l. We followed Kochanek|( |2004| l in setting the present- 
day velocity dispersion of galaxies cr pec to 235 km s _1 . The 
function f(z) is the cosmological growth factor; we approxi- 
mate it as / ex J7 M (z) 0,6 . The stellar velocity dispersion a L of 
the lens we estimate from the monopole compone nt b of our 
lens m odel (see Section E\ using the relation (e.g., Kochanek 
[2006l > 



b = 4iT 



c 2 £>os 



(4) 



All velocities are corrected for cosmological time dilation. 
The factors of two convert one-dimensional velocity disper- 
sions to two dimensions. Multiplying the estimated source 
speed vj_ by the delay between our optical/IR observations 
and their X-ray counterparts, we determined the distance trav- 
eled by the source along the caustic pattern, and from the 
structure function read off the predicted error in our X-ray 
flux ratio. 

For each X-ray flux ratio, we added in quadrature the un- 
certainty contributions from measurement errors (including 
extra uncertainty for image A2 of WFI J2026 because of its 
degeneracy with the PSF width), intrinsic quasar variability, 
and microlensing variability. These uncertainties are given in 
Table 15] In general, the largest contributions were systematic 
uncertainties from microlensing variability, though the statis- 
tical measurement errors were also substantial in some cases. 
As in the optical case, the uncertainty due to quasar variability 
was relatively insignificant. 
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Table 6 

Lens Model Parameters 





Primary Lens 


Secondary Lens 




Magnification 3 




Quasar 


b 


7 


07 " 


b 2 


*2 C 


V2 C 


HM 


HS 


LM 


LS 


HE 0230 


0."87 


0.11 


-60.°0 


0"33 


-0.283 


+0.974 


+9.42 


-9.65 


+4.95 


-1.35 


MGJ0414 


l."14 


0.11 


+74°6 


0."12 


-0.385 


+1.457 


+22.9 


-24.2 


+6.23 


-3.11 


HE 0435 


1"20 


0.078 


-13.°8 








+7.49 


-7.90 


+7.14 


-4.73 


RXJ0911 


0"97 


0.27 


+7.°2 


0."24 


-0.754 


+0.665 


+11.0 


-5.96 


+1.97 


-4.99 


SDSS J0924 


0"87 


0.063 


+84.° 8 








+14.9 


-13.0 


+6.62 


-6.55 


HE1113 


0"33 


0.040 


+37.°7 








+15.8 


-16.7 


+12.6 


-9.59 


PG1115 


1"03 






2"57 


-10.866 


-5.300 


+19.7 


-18.9 


+5.09 


-3.37 


RXJ1131 


1"86 


0.16 


-73.° 6 








+13.2 


-22.7 


+12.6 


-1.05 


SDSS Jl 138 


0"67 


0.10 


+32°6 








+7.17 


-6.68 


+5.17 


-3.64 


SDSSJ1330 


0"94 


0.16 d 


-32.° 2 d 








+27.1 


-27.2 


+8.41 


-5.50 


WFIJ2026 


0"66 


0.11 


-90.°0 








+13.7 


-11.5 


+3.78 


-4.01 


WFIJ2033 


V.'Ol 


0.11 


+36.°0 


0."25 


+0.229 


+2.020 


+6.04 


-3.80 


+3.88 


-2.46 



[1 HS: highly magnified saddle point; HM: highly magnified minimum; LS: less magnified saddle point; 

LM: less magnified minimum. 

b Position angle of external shear 7. All angles are measured in degrees east of north. 

c Fixed position of secondary galaxy, relative to main lensing galaxy, in arcseconds. Allowed to vary 

radially in the case of PG 1115; see the text. The positive directions of x and y are west and north, 

respectively. 

d Ellipticity and position angle of the lens galaxy for the SIE model of SDSS J 1330. 



5. LENS MODELS 



We used the lensmodel program of Keeton (2001 ) to create 
parametric models of each lens. The models were constrained 
by the positions of the four lensed images, and that of the 
lensing galaxy, a total of 10 constraints. We did not use fluxes 
for constraints, since most of our lenses suffer from flux ratio 
anomalies. Nor did we use time delays. The positions for the 
quasar images and lens galaxies came from the CASTLES 
survey. 

Our default model consisted of a singular isothermal sphere 
for the lensing galaxy, with a quadrupole component of the 
potential provided by a constant external shear. With the po- 
sition of the lens fixed, this model had five free parameters: 
the monopole strength of the lens, the magnitude and direc- 
tion of the shear, and the position of the source. 

In some cases, such as HE 0435, this simple model fit the 
image positions well. But in several cases the x 2 goodness 
of fit was poor enough to warrant further complexity in the 
model. In these cases, we made changes to the model moti- 
vated by the appearance of the lens galaxy. 

The lens galaxy of HE 0230 has a prominent companion, 
located close to image D. We modeled this companion as a 
second isothermal sphere, fixing its position to its HST mea- 
sured value but allowing its mass to vary. Despite only adding 
one free parameter to the model, this addition improved the fit 
considerably. We followed a similar strategy for MGJ0414, 
RXJ0911, and WFIJ2033, each of which displays a faint 
smudge in its HST image which is arguably a satellite to the 
lens galaxy. Adding secondary lenses at the positions of the 
smudges improved the goodness of fit to acceptable levels. 

The lens galaxy in PG 1115 does not have a nearby compan- 
ion, but is a member of a small galaxy group cent ered to the 
southwest o f the lens. Following the example of |Schechter| 
|et al.| (JT997J, we explicitly modeled the group as a second 
isothermal sphere for this lens. We parameterized its position 
using polar coordinates, and allowed its mass and distance 
from the main galaxy to vary while fixing its position angle to 
that of the brightest galaxy in the group. We did not include 
an external shear in this fit. 



Finally, the lens galaxy in SDSS J1330 displays significant 
ellipticity. In this case we used an isothermal ellipsoid instead 
of including an external shear. We allowed the ellipticity and 
position angle to vary, along with the galaxy position, for a 
total of seven free parameters. Since there are no HST data 
for this lens, we used our measured image positions for con- 
straints (see Section[3]and Table[3]l. SDSS J 1330 was the only 
lens system where an isothermal ellipsoid made for a better fit 
than an isothermal sphere with external shear. 

The salient features of our best-fit models are listed in Ta- 
ble [6] We used the components of the magnification tensor 
predicted by these models (namely, the local convergence k 
and shear 7) to generate th e mi crolensing magnification pat- 

We also used the predicted 



terns described in Section 6.1 



scalar magnifications to estimate the microlensing-free flux 
ratios, and assumed that any discrepancies between the pre- 
dicted and observed flux ratios was due to microlensing. We 
believe that this is a good assumption for our sample of lenses 
because there is not strong evidence for other potential ex- 
planations. For example, differential extinction is unlikely 
to be a fact or in these early-type lens galaxies ( Kochanek & 
Dalal|[2004] l; this is supported b y the general lack o f strong 
hydrogen absorption in X-rays ( |Pooley et al.||2007) . Like- 
wise, the chromatic variations we see in the flux ratios ar- 
gues against millilensing by dark matter substructure. In the 
case of MGJ0414, however, mid-infrared flux ratios mea- 



sured by Minezaki et al. (2009) disagree with the predicted 
ratios. Since the dusty torus region where this light originates 
is too large to be affected by microlensing, it is likely that 
millilensing is affecting this lens. In this case, we adopted the 
reported mid-IR flux ratios as the "baseline" rather than the 
model ratios, attributing any optical departures from them to 
microlensing. 

The predicted flux ratios are model-dependent to some ex- 
tent; this systematic error is difficult to quantify, but ought 
to be sma l l com pared to the other uncertaint ies. |Dalal & 
Kochanek (2002) suggest error bars of 10%; Keeton et al. 
( 2003 ) call this "quite conservative." We estimated the model 
uncertainties in the flux ratios as 0.05 mag for well-separated 
quasar images, or 0.03 mag for members of a close pair 
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or triplet, since these configurations have model-independent 
asymptotic flux ratios ( |Keeton et al.||2003 2005 and refer- 
ences therein). We added this uncertainty in quadrature to 
those listed in Table [5] 

6. BAYESIAN SOURCE SIZE ESTIMATION 

The effect of a finite source size is to smooth out the tem- 
poral variations caused by microlensing. The movement of 
a large source across a network of microlensing magnifica- 
tion caustics is equivalent to a point source moving across a 
version of the caustic map that has been smoothed by convo- 
lution with the surface brightness profile of the source. Since 
microlensing variability is stronger in X-rays than at optical 
wavelengths, it is thought that the X-rays come from a very 
compact region, whereas the optical source is much larger 



dPooley et al.|2007l[2009)|Morgan et al.|2008||Dai et al.|2010| >. 
For single-epoch observations like ours, microlensing magni- 
fications are more naturally treated as flux ratio anomalies, 
where the flux ratios of the quasar images fail to match those 
predicted by lens models, than as variability. We expect the 
flux ratio anomalies to be stronger in X-rays than at optical 
wavelengths, and stronger in blue light than in red light, most 
of the time. However, different regions of the magnification 
maps have different responses to a change in the source size. 
Indeed, in some regions (e.g., at the center of an astroid caus- 
tic), convo lution by a lar ge source profile will increase the 
anomaly. |Poindexter et al. ( 2008} saw an example of this in 
the doubly lensed quasar HE 1 104- 1805. 

Single-epoch measurements are much less expensive than 
long-term monitoring campaigns, especially with multiple 
wavelengths. But they have the disadvantage that they are un- 
able to constrain the local characteristics of the microlensing 
magnification map, leading to a degeneracy between the natu- 
ral range of possible magnifications and source size effects. A 
strong anomaly is evidence for a relatively small source, but a 
weak or absent anomaly could be due to a large source or to 
the chance location of a small source on a relatively calm por- 
tion of the magnification map. PG 1 1 15 is a real-life example 
of this: for several years its optical A2/A1 ratio maintained 
a low-level anomaly that was revealed to be due to a large 
source size when the same ratio was found to be extremely 
low i n X-ra ys (Pooley et al. 2006). Our analysis method fol- 
lows |Pooley et al.| ( ]2007| l in using X-ray flux ratios to break 
this degeneracy, but in a more quantitative way. 

A similar uncertainty exists when using size estimates at a 
range of wavelengths to estimate the wavelength dependence 
of the size of quasar accretion disks. Because of the un- 
known response of the microlensing magnification to changes 
in source size, a single-epoch observation of a single lensed 
quasar cannot tell us much about this dependence. We stress, 
however, that it does allow for estimates with meaningful er- 
ror bars. The way to decrease this uncertainty is to sample a 
greater fraction of the caustic map, either via time-series ob- 
servations of a single lens or by observing a sample of lenses. 
We have adopted the latter approach. 

We have developed a quantitative Bayesian method for esti- 
mating the angular size of quasar accretion disks at a range of 
optical and IR wavelengths using our measured optical/IR and 
X-ray flux ratios. This approach relates the posterior probabil- 
ity distribution for the source size to the likelihood of observ- 
ing the reported X-ray flux ratios, assumed to originate in a 
very compact source, simultaneously with the observed opti- 
cal/IR ratios from a larger source (see Equation (|5]l). The bulk 



of this section is devoted to the calculation of these likelihood 
distributions. 

6.1. Microlensing Magnification Patterns 

We simulated the microlensing of a finite-size source using 
magnification p atterns created u sing the inverse ray-tracing 
software of Wambsganss et al. (1990i. Each pattern corre- 
sponds to a specific quasar image, and was constructed using 
the local convergence n tol and shear 7 specified by our lens 
models (see Section Bl). The convergence (which is propor- 
tional to the surface mass density) we divided into k c due to 
continuously distributed matter (e.g., dark matter) and k* due 
to stellar-mass point lenses. For each quasar image we cre- 
ated 12 patterns, varying the surface mass fraction in stars 
k* = «*/«tot from 1.47% to 100%, in logarithmic steps. 

The point mass microlenses were scattered randomly across 
the microlensing patterns, with masses drawn at random from 
a fixed mass function. The primary effect of the mass func- 
tion chosen for the simulated microlenses is to set the physi- 
cal size of the microlens Einstein radius, the natural scale for 
microlensing simulations. Different mass functions also lead 
to differently shaped magnification histograms, but this effect 
is less important. We were not able to make an estimate of 
the mean mass (m) using our single-epoch measurements, so 
we adopted a fixed broken power-law m ass function nearly 
identical to the well-known |Kroup"a ( 2001) > initial mass func- 
tion. Between 0.08M Q and 0.5M Q , its logarithmic slope was 
—1.8; above 0.5M Q it steepened to -2.7. We cut off the mass 
function at 1.5 Mq, because the stellar populations in these 
early-type lens galaxies are typically old. With this mass func- 
tion, the average microlens mass (m) was 0.247M Q . With this 
fixed mean mass came a fixed transformation between the an- 
gular sizes we measure (in units of Einstein radii) and phys- 
ical sizes at the distance of the quasar (in units of cm); any 
adjustment to it will scale our measured disk sizes according 
to(m) 1 / 2 . 

Each pattern encodes the deviation in the magnification of 
its quasar image from that produced by a smooth mass dis- 
tribution, as a function of the position of the source. They 
are 2000 pixels on a side, and represent a square region of 
the source plane with a side length 20 times the projection 
of the Einstein radius of a solar-mass star in the lens galaxy. 
This distance is roughly 5 x 10 17 cm for the lenses in our sam- 
ple, though the exact numbers depend weakly on the lens and 
source redshifts. With these dimensions, the patterns are large 
enough that their histograms are insensitive to the exact con- 
figuration of the stars. The pixel size is 2.5 x 10 14 cm, or a 
few gravitational radii for a 10 9 M Q black hole. This is much 
smaller than the size of the optical a ccretion disk a nd is com- 
parable to the source size in X-rays (Morgan et al. 2008 ; Dai 
|etal.|2010) . 

We modeled sources of finite size by convolving the mag- 
nification patterns with circular Gaussian kernels of varying 
half-light radii. Though a Gaussian is not a very physical 
choice for a source profile, it has been shown (Mortonson 
|et~aT1|2005| |Congdon et aLl|2007> that the half-light radius 
is much more important when simulating finite sources than 
are the details of the radial profile. This convolution blurs the 
patterns and reduces their dynamic range, causing their his- 
tograms to become narrower. The result of our efforts was 
a family of patterns for each quasar image, parameterized by 
stellar mass fraction and source size. 
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The local surface mass fraction in stars as a function of 



position within a galaxy is poorly constrained. Pooley et al. 
j2009| ) used X-ray flux ratios of PG 1 1 15 to estimate the stel- 
lar mass fraction at the locations where the quasar images in- 
tersect the lens galaxy. Because the quasar appears to be com- 
pact in X-rays (Pooley et al.|2007[|Chartas et al72009), they 



were able to safely neglect the finite source size, and found 
a probability distribution for the stellar mass fraction k* that 
was peaked near 10%. Since a detailed study of stellar mass 
fractions is beyond the scope of this work, we marginalized 
over th is parameter usin g t he re sults found by Pooley et al. 
( |200"9l >(see Sections [O] and |63]>. 

6.2. Size Estimation Using Optical Ratios 

We used a simple Bayesian method to calculate the pos- 
terior probability distribution for the half-light radius of the 
source ri i 2 at each wavelength. Assuming a uniform prior for 
the source size, the probability is 



1/ 2 |{wj}) oc£({m,}|r I/2 ) 



(5) 



where {m,} is the set of quasar image magnitudes, i runs from 
1 to N[ m = 4, and C is the likelihood of having observed the 
magnitudes {m,-} given a source half-light radius. Each mag- 
nitude mi may be expressed as the sum of the unmagnified 
quasar magnitude ms, macro- and micro-magnification terms 
/!,■ and A/i,, and measurement noise n,\ 



nil = tns + A*; + A/i; + n,- 



(6) 



The source magnitude nis is not observable, and because we 
did not flux-calibrate our measurements, it is in uncalibrated 
units. Since we ultimately marginalized over m s , these limi- 
tations are unimportant. 

We computed the right-hand side of Equation (T5]l by com- 
paring the measured image fluxes mi to histo grams of the 
magnification patterns described in Section |6Tj Suitably nor- 
malized, each histogram was an estimate of the distribution of 
microlensing magnifications: 



P;(A/^|r 1/2 ,fc*) 



(7) 



where the pattern sourcing the histogram was created using 
the stellar mass fraction k* and was convolved with a source 
with radius r\ji. 

Figure [2] shows an example of four such magnification his- 
tograms, one for each image in SDSS J 1330, for a particular 
source size and stellar mass fraction. The abscissa of this plot 
denotes the microlensing deviation A/i in magnitudes. Atop 
the histograms have been plotted vertical lines, indicating the 
observed deviation from the model in the if band. In the case 
of no microlensing, the lines would be co-located at zero. The 
presence of four vertical lines misleadingly implies more in- 
formation than is available; in fact, we have not measured the 
four image magnifications, but only their three ratios. Thus, 
the average position of the four lines is unknown. In Figure[2j 
we have set it to zero, but in reality the four lines are free to 
slide from side to side in formation. Each horizontal posi- 
tion corresponds to a distinct value of the quasar's unknown 
magnitude ms- 

We shifted each histogram horizontally to obtain /?,( A/z, - 
fWj + Mi) = Pi(~ m s~ n d- These shifts brought the four vertical 
lines in Figure [2] together. Because the shifts depended on 
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Microlensing magnification (magnitudes) 

Figure 2. Histograms derived from magnification patterns for the four im- 
ages A, B, C, and D of SDSS J1330 (top to bottom), with a source radius 
3.2% of a solar-mass Einstein radius and a stellar fraction of 32%. Flux ra- 
tios from our i' band data have been overplotted in dotted lines. The abscissa 
denotes deviation from the predicted flux ratios; if the dotted lines all fell at 
zero, the lens would have no flux ratio anomaly. Magnifications are in mag- 
nitudes, so positive numbers denote demagnification. The histograms of the 
negative parity images (B and D) are notably wider than those of the positive 
parity images. The same is true of the more highly magnified images (A and 
B). 

the measured fluxes, we marginalized over the measurement 
errors: 



Pi(-m s ) = 



2ira; 



dn, exp 



2a} 



pi(-m s -ni) . (8) 



This amounted to a convolution of the shifted histograms with 
the normal distributions implied by the measurement uncer- 
tainties er,. These shifted and convolved histograms were 
probability density functions of the quasar magnitude nis- We 
multiplied them together and integrated their product to find 
the likelihood of measuring all the observed flux ratios, inde- 
pendent of the quasar flux: 



£({m,-}|ri/2,&*) = / dm S (j)(m s ) 



Y\_Pi(-ms) 



Li=l 



(9) 



where N\ m = 4 and 4>(nis) is the prior on the luminosity of the 
quasar. For this luminosity function we used a power law of 
slope -2.95, as estimated for bright quasars by Richards et al. 
(2006). Our use of a scale-free luminosity function obviated 
any difficulties due to the uncalibrated quasar flux. 

Finally, we marginalized over the stellar mass fraction k* to 
find the right-hand side of Equation Q: 



^{{ m i}\ r \/i) = j dk*p(k*)£({mi}\r l/2 , 



(10) 



For our prior p(k*) on the stellar mass fraction we interpo- 
lated the estimat ed probab ility distribution for its value given 
in Figure 6 of Pooley et al. ( 2009[ ). It should be noted that al- 
though we have written the marginalizations in Equations §8\- 
( [T0| as integrals for clarity, they were actually calculated as 
discrete sums. 
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log(ri/ 2 /rEin) 

Figure 3. Posterior probability distribution for the size of SDSS J 1 330, using 
the i' band data. The histograms were convolved with the measurement un- 
certainties, and the quasar luminosity function was used as a weighting factor 
when integrating the product of the histograms. A logarithmic prior was used 
for the source size. 



We calculated this relative likelihood for 3 1 half-light radii, 
spaced logarithmically between 0.01rE; n and lOre;,,, where 
rEin is the Einstein radius of a solar-mass microlens. We let the 
posterior probability of the source size equal this likelihood, 
implicitly adopting a logarithmic prior on the source size (be- 
tween these limits; zero elsewhere). For completeness, we 
also used a linear prior; this simply involved multiplying the 
likelihood distribution by a power law of unit slope. 

Figure[3]shows the posterior distribution resulting from the 
use of this technique on the /'-band flux ratios of SDSS J 1330. 
This figure demonstrates the ability of this method to rule 
out large source sizes (e.g., greater than an Einstein radius). 
The decrease in probability to zero occurs when the grow- 
ing source size causes two of the shifted histograms to be- 
come narrow enough that they no longer overlap at all; at this 
point the likelihood is zero. Measurement errors widen the 
histograms via the convolution in Equation dS), causing the 
probability distribution to widen as well. There is no corre- 
sponding decrease in probability for small source radii. This 
is to be expected, since the optical flux ratio anomalies sim- 
ply indicate that the source is undergoing microlensing. To set 
a lower limit on the size, we must discern why the observed 
anomalies are not stronger than they are — is the caustic map 
being smoothed out by a large source, or is the source simply 
lying in a non-perturbed region? Though some size discrim- 
ination is possible by comparing the posterior distributions 
derived from blue filters versus red filters, a more effective 
method to rule out very small sizes is to compare the optical 
flux ratios to those in X-rays. 

6.3. Size Estimation Using Optical and X-ray Ratios 

SDSS J1330, the lens we chose for the examples in the pre- 
vious section, is unique among our sample in that it has never 
been observed in X-rays. But for the remaining lenses in our 
sample there is at least one measurement of the X-ray flux ra- 
tios. Since the microlensing variability is in general stronger 
in X-rays than at optical waveleng ths, the X-rays must come 
from a quit e compact region (Pooley et al. 2007 2009 Dai 
et al. |2Q10) >. Any decrease in the anomalies in the optical 
and IR must be due to the finite size of the source at these 



wavelengths. Comparing the optical flux ratios to X-ray ra- 
tios gives us the ability to set a lower limit on the size of the 
optical-emitting accretion disk. 

We again wished to find the value of the right-hand side of 
Equation d3), but this time, we wanted the likelihood of ob- 
serving both the X-ray and optical flux ratios simultaneously, 
given a source size. Because the X-rays are thought to origi- 
nate from a very compact region, we assumed that the X-ray 
ratios were drawn from the original (unconvolved) magnifica- 
tion pattern. Then, for each image, we constructed a family 
of two-dimensional histograms: 



Pi (At4,An°\r 1/2 ,k* 



(11) 



where A/i? and A/if indicate values taken from the con- 
volved (optical) and unconvolved (X-ray) versions of the 
same pattern. These histograms, after appropriate scaling, are 
joint probability distributions for the optical and X-ray mi- 
crolensing magnifications; they are the two-dimensional gen- 
eralization of the histogram in Equation dTli. Figure HI shows 
two-dimensional histograms corresponding to the images of 
PG 1115, with the X-ray magnifications drawn from the orig- 
inal pattern and the optical magnifications drawn from the 
same pattern smoothed by a source of varying half-light radii. 
Histograms for all four images are shown for a single source 
size (0.01r E ; n ), and the histogram for the highly magnified 
saddle point image A2 is shown for a variety of source sizes. 
The maps used for these histograms all had a stellar mass frac- 
tion k* =0.1. 

We shifted the histograms along the optical and X-ray axes 
in a manner analogous to the one-dimensional case, using 
the measured flux ratios in X-rays and in the chosen opti- 
cal/IR filter. For very small source sizes, the two patterns 
from which each histogram was constructed were nearly iden- 
tical, so the histograms showed a significant correlation be- 
tween A/i° and A/z x . Thus, for a very small source size to 
be assigned any significant probability, the optical and X-ray 
shifts had to be the same. As the optical source size increased, 
the histograms' extent in the A/i° direction decreased, as did 
the correlation between the magnification in the two bands, 
until for very large source sizes the histograms were nearly 
aligned with the A/i x axis. In this case, the optical shifts had 
to be nearly zero for the probability to be appreciable (see 
Figure HI. 

In the one-dimensional case, we propagated the uncertain- 
ties in the measured flux ratios by convolving our shifted 
histograms with Gaussians representing the uncertainties. 
In this two-dimensional context, convolution (e.g., with a 
two-dimensional Gaussian) gave us insufficient information. 
Since the source size likelihood in every filter depended on 
the same X-ray flux ratios, the X-ray uncertainties led to cor- 
related errors in the size estimates at different wavelengths. 
So we instead performed a Monte Carlo integration in order 
to propagate the full covariance matrix through the analysis 
described in the next paragraph. We describe the error propa- 
gation in the paragraph that follows. 

We multiplied the shifted histograms and marginalized over 
the unknown X-ray and optical source fluxes to find the like- 
lihood of simultaneously measuring the observed X-ray and 
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-4-2 2 4-4-2 2 4 
Optical magnification (magnitudes) 

Figure 4. 1<t, 2a, and 3cr contours of the two-dimensional microlensing 
magnification histograms for the four images of PG 1115, after the magnifi- 
cation map has been convolved with Gaussian source profiles with varying 
optical half light radii. In the upper left panel, histograms for all four images 
are shown, with an optical source size of 10~ 2 rEi n - The other panels show the 
evolution of the histogram of image A2 as the optical size increases; the half 
light radii are labeled in the upper left of each panel, in units of r^m ■ As the 
source gets larger, the correlation between X-ray and optical magnifications 
decreases, and the histogram narrows in the optical direction. For even larger 
source sizes, the histograms get very narrow. 

optical fluxes: 



£({mf ,m°}\r 1/2 ,k*) = / dm x s dm°(, 



n 



Pi(-m^,-m^) 



(12) 



In this case, the weighting function 0(mf , m°) incorporated a 
prior both on the luminosity function of quasars and on the 
correlation between their X-ray and optical fluxes. The latter 
is often parameterized using aox = 0.3838 log^kev/^sooA)- 
Gibson et al. (2008) find that aox = _ 0-2171og(L 2 5ooA) Pl us 
a constant olrset, with a scatter of 0.1. This implies that 
m* = 0.4346ra5? plus a constant offset, with a scatter of 0.65 
mag. So we constructed a band centered on a line with a 
slope of 0.4346 in the (mg,m?) plane, and with a Gaussian 
cross-section. In order to be conservative, we doubled the ob- 
served scatter and set the width of the Gaussian to correspond 
to a = 1 .3 magnitudes. We adjusted the constant offset so that 
the line passed through the mean of the product histogram. 
We multiplied this function by a power law of slope -2.95 (as 
in the one-dimensional case) in the direction along the line of 
correlation to represent the luminosity function. We note that 



1.5 
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Figure 5. Posterior probability distribution for the size of PG 1115 in the 
i' band, resulting from considering both / '-ban d and X-ray flux ratios. The 
Monte Carlo method described in Section l6~3l was used to account for mea- 
surement uncertainties. A logarithmic prior was used for the source size. 



the slope of the luminosity functio n of X-ray sele cted bright 
quasars is very similar at -2.8 (Silverm an et al.|2008) >. As be- 
fore, we repeated this technique using magnification patterns 
representing a range of stellar mass fractions, and marginal- 
ized over this parameter using a weighted average of the re- 
sulting likelihood distributions (see Equation ( 10 1). 

The Monte Carlo integration which allowed us to propa- 
gate the measurement uncertainties into the likelihood distri- 
butions consisted of repeating the process of shifting, multi- 
plying, and integrating under the histograms for 1000 sets of 
flux ratios drawn from Gaussian distributions implied by the 
measurement uncertainties. This gave us 1000 distributions 
for the disk size in every filter. From these we took the most 
likely values (the modes) and constructed the covariance ma- 
trix Cij = Cov(r 1 /2,;,'"i/2j) 5 where i and j refer to optical/IR 
filters. The diagonal elements of the matrix are simply the 
variances of the half-light radii due to measurement uncertain- 
ties. To this matrix we added a second, diagonal matrix 5y, 
which expressed the "intrinsic" width of the first of our 1000 
likelihood distributions, the one derived from the reported flux 
ratios. This intrinsic width arose not from measurement un- 
certainties, but simply from the widths of the magnification 
histograms. The diagonal elements of S were the squares 
of the halved 68% confidence intervals of these distributions, 
and the remaining elements were zero. We describe in the 
next section how we used the inverse of the sum of these co- 
variance matrices 



W = (C+ST 



(13) 



as a weight matrix in a chi-square fit. 

The fact that the optical flux ratios demonstrate microlens- 
ing magnifications allows us to place upper limits on the size 
of the optical disk, since a very large disk will be unaffected 
by microlensing. Likewise, the fact that the X-ray flux ratios 
are different from (typically, larger than) the optical ratios al- 
lows us to place lower limits. For small optical sizes, the dif- 
fering flux ratios cause the highly correlated histograms to be 
shifted in a direction different from that of the correlation, 
lowering the likelihood (see Figure HI. The posterior proba- 
bility distribution for the /'-band radius of the accretion disk 
of PG 1115, determined in this manner, is shown in Figure 
and shows a central peak, with a decrease in probability to 
left due to the interaction of the X-ray and optical flux ratios, 
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Rest Wavelength (microns) 

Figure 6. Medians (squares) and modes (crosses) of the probability distri- 
butions for the half-light source radius. The error bars correspond to the 
quadratic sum of the intrinsic width of the distributions and the scatter from 
the Monte Carlo propagation of measurement uncertainties. Black (gray) 
points indicate a logarithmic (linear) prior. The solid line is the prediction for 
the half-light radius (in cm) of the standard thin disk model, for the estimated 
black hole mass. 



and a decrease to the right due primarily to the optical flux 
ratios alone. 

We tested the accuracy of our analysis method by apply- 
ing it to artificial flux ratios generated using the predictions 
of the thin disk model. We calculated the half-light radii that 
this model predicts for PG 1 1 15 in each filter, then convolved 
the magnification patterns appropriate for the four images of 
this lens (with stellar mass fraction k* = 0.1) with source pro- 
files of these radii. Then we drew microlensing magnifica- 
tions from 30 random locations in these patterns. We drew 
X-ray ratios from the same locations in the corresponding 
non-smoothed patterns. We ran our analysis on these 30 sets 
of simulated flux ratios, with uncertainties equal to those as- 
signed to the real PG 1115 data, to check if we recovered the 
input temperature profile. We fit a line to the recovered sizes 
to determine the logarithmic slope and normalization of the 
wavelength-dependent half-light radius, using the procedure 
described in Section 17] On average, the best-fit overall nor- 
malization log(ri/2/cm) was 16.1, with a scatter of 0.2 dex, 
compared to the input value of 15.7. The average output slope 
was 1.2, with a scatter of 0.1, versus the input value of 1.33. 
The general agreement of the recovered parameter values with 
the input values indicates that the analysis method produces 
sensible results. The small offset between the input and out- 
put normalizations indicates that our analysis may overesti- 
mate the sizes slightly. 



7. RESULTS 

The result of the analysis described in Section[6] was a pos- 
terior probability distribution for the half-light radius of the 
quasar accretion disk in each of our 12 lensed quasars, in each 
filter. The medians and modes of these distributions are listed, 
along with lc confidence intervals, in Table [JJ and plotted in 
Figures [6]-[9] as a function of rest wavelength. The error bars 
are taken from the diagonal elements of the covariance ma- 
trix W _1 defined in Equation ( |13| ). The choice of prior for 
the half-light radii has a smalleffect on the measured val- 
ues. In general, the linear prior favors values of r\n ~ 0.2 dex 
larger than the logarithmic prior does. This is to be expected, 
since the linear prior has more probability density at larger 
sizes. The slope with wavelength is generally unaffected by 
the choice of prior. Since a logarithmic prior is more appro- 
priate for a scale-free quantity such as the half-light radius, 
we adopt these values. 

The half-light radii are measured in terms of the Einstein ra- 
dius of a solar-mass microlens r-£ m - We converted these angu- 
lar sizes to physical radii in the source plane using the angular 
diameter distances Dos, Do L , and D L s (observer-to-source, 
observer-to-lens, and lens-to-source, respectively). The phys- 
ical source radius depends weakly on the redshifts of the 
source and lens, and the chosen cosmology. It also depends 
on our choice of the stellar mass fun ction , going as the square 
root of the mean mass (see Sectio n |6.1[ l. We show physical 
distances on the left axis of Figures lofl^J with the correspond- 
ing fraction of an Einstein radius on the right axis. 

In order to compare the measured source sizes to the thin 
disk predictions, we must estimate the Eddington fraction, 
accretion efficiency, and black ho l e mas ses of our sample of 
quasars. Following |Po oley et al.| (|20 07|), we set fay = 0.25 

77 = 0.i5 dYu_ 



( |Kollmeier et al.|2006[ )"and~ 



QYu & Trerna ine 2002). 
We turned to the literature to find black hole mass estimates 
for most of our sample. For MG J0414, HE 0435, RX J091 1, 
SDSSJ0924, PG 1115, andR XJ1131 we adopted the virial 
mass estimates of |Peng et al| (|2006), which are based on the 
widths of the broad emission lines Civ, MgH, and H/3. For 



SDSS Jl 138 we adopted that of M organ et aLlpHO) !, calcu- 
lated using the same m ethod. We adopted the masses esti- 
mated by Pooley et all l |2007) for HE 0230 and WFIJ2033 
using their optical and X-ray luminosities. For the remaining 
quasars (HE 1113, SP SS J1330, and WFI J2026) we used the 
luminosity method of Pooley et al. (2007), inferring a bolo- 
metric luminosity from the optical flux. Trie 
used the relation 



le optical estimates 



Lbo\ = 9[Xfx] 51QoA 4Trdf 



(14) 



(Kaspi et al|2000| l, where d[ is the luminosity distance to the 
quasar, and the rest-frame 5 100 A flux is extrapolated from the 
magnification-corrected HST NICMOS F160W flux of the 
LM image using an a ssumed power-law spectrum f\ ~ A -1,7 
( Kollmeier et al|2006) . We used our assumed Eddington ra- 
tio /Edd = 0.25 to calculate black hole masses. The luminosity 
and virial mass estimates for our sample of quasars are listed 
in Table [8] We adopted the virial estimates where available, 
and the luminosity masses elsewhere. We note that in cases 
where both are available, the luminosity-based estimates are 
systematically smaller than the virial estimates by ~ 0.5 dex; 
this could lead to under-predictions of the accretion disk size 
by factors up to ~ 2 for the cases where we use the luminosity 
estimate in the absence of a virial estimate. This systematic 
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Figure 7. Medians (squares) and modes (crosses) of the probability distribu- 
tions for the half-light source radius. The error bars correspond to the sum in 
quadratures of the intrinsic width of the distributions and the scatter from the 
Monte Carlo propagation of measurement uncertainties. Black (gray) points 
indicate a logarithmic (linear) prior. The solid line is the prediction for the 
half-light radius (in cm) of the standard thin disk model, for the estimated 
black hole mass. 
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Figure 8. Medians (squares) and modes (crosses) of the probability distribu- 
tions for the half-light source radius. The error bars correspond to the sum in 
quadratures of the intrinsic width of the distributions and the scatter from the 
Monte Carlo propagation of measurement uncertainties. Black (gray) points 
indicate a logarithmic (linear) prior. The solid line is the prediction for the 
half-light radius (in cm) of the standard thin disk model, for the estimated 
black hole mass. 



bias is a little surprising, since the bolometric luminosity tech- 
nique is calibrated using the virial method (Kollmeier et al. 
2006). It may be partially due to a bias in our magnification 
correction, but it is difficult to imagine that we have system- 
atically overestimated magnifications by a factor of three. 

Using these black hole mass estimates, we plot the predic- 
tion of the thin disk model for each quasar in Figures [6}j9] 
Comparing the microlensing sizes to these curves, it is im- 
mediately apparent that the measured sizes are systematically 
larger than predicted by the thin disk model and vary more 
weakly with wavelength than predicted. 

For each lens we fit a power-law model Y(X) to the medians 
of the distributions: 



Y(\) = Y c + vlog 



(15) 



where the central wavelength A c is the geometric mean of the 
rest wavelengths and Y c is the logarithm of the half-light ra- 
dius of the source at that wavelength, in cm. Since the errors 
on the medians are correlated, we minimized the statistic 



X 



= (y-Y(X)) J -W-(y-Y(X)) 



(16) 



where y is the vector of median sizes measured at wavelengths 
A, and W is the weight matrix defined in Equation ( 13 I, The 
parameters Y c and v of the best-fit power laws are listed in 
Table 



In nearly every case, the measured disk size at the central 
wavelength \ is larger than the thin disk model predicts, by 
factors ranging from ~ 2 to more than 30. The logarithm of 
this factor is listed for each quasar in Table [9] Comparing 
the logarithmic offset to our errors in Y c , we have ruled out 
the thin disk normalization by at least 3cr in all cases except 
MG J0414, and as ma ny as 10cr in s ome cases. This result is 
consistent with that of Pooley et al. (2007), though our anal- 
ysis method is more quantitative. The average offset between 
Y c and the thin disk prediction at the same wavelength is 0.89 
dex, or a factor of 7.5. A chi-square test comparing the mea- 
sured sizes to predicted sizes yields xV-^W = 46. 

In Figure [lOj we plot the half-light disk radii of our sample 
(excluding SDSS J1330) at fixed rest wavelength versus the 
black hole mass. The disk radii are calculated using the best- 



fit model (see Equation ( 15 i) at A = 25 00 A. This wavelength 
is chosen to match the one used by Morg an et al.| ( |2010[ ), 
who extrapolate from single-wavelength microlensing obser- 
vations to create a similar plot. The errors on the half-light 
radii are the uncertainties on Y c (see Table 9l, and uncertain- 
ties of 0.4 dex have been assigned to the black hole masses 
(but not plotted). We multiplied the mass uncertainties by the 
thin disk model slope of 2/3 and added them in quadrature to 
the uncertainties on the half-light radii. With these uncertain- 
ties, the best-fit slope for r\n versus Mbh is 0.27 ±0.17. The 
X 2 for this fit is 17.3 for 9 degrees of freedom. Thus, we find 
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Table 8 








Black Hole Mass Estimates 






^bol.opt 


ibol.X 


MBH.opi 


MBH.vir 


Quasar 


(10 46 erg s" 


) (10 46 erg s" 


) (1O 9 M ) 


(10 9 M o ) 


HE 0230 


0.29 


0.63 


0.092 




MGJ0414 


3.6 


2.8 


1.1 


1.82 


HE 0435 


0.38 a 


0.46 a 


0.12 


0.50 


RXJ0911 


1.3 


1.3 


0.41 


0.80 


SDSSJ0924 


0.06 


0.03 


0.019 


0.11 


HE1113 


0.27 a 


0.10 a 


0.087 




PG1115 


1.1 


0.66 


0.35 


0.92/1. 23 b 


RXJI13I 


0.08 


0.13 


0.025 


0.06 


SDSSJ1138 


0.38 a 


0.25 a 


0.12 


0.04 c 


SDSSJ1330 


4.7 a 




1.5 




WFIJ2026 


2.5 a 


1.1" 


0.79 




WFIJ2033 


0.57 


0.38 


0.18 




Note. — 


Unless other 


wise indicated 


bolometric 


luminosity 


estimates are 


from Pooley 


et al. (2007 1 and virial mass estimates 



are from Peng et al. (2006). 

a This work. 

b Two values are from the C IV and Mg II lines, respectively. We 

adopt the Mg II value. 

c |MorganetaLl ( (20T0) 



Table 9 

Best-fit Parameters for Disk Size Versus Wavelength 





A c a (A) 




Logarithmic Prior 






Linear Prior 




Quasar 


Y c 


logf/r^Apred) 


V 


Y c 


logOl^Apred) 


V 


HE 0230 


2763 


16.57±0.15 


+1.71 


-0.56 ±0.47 


16.79±0.18 


+1.93 


-0.36 ±0.43 


MGJ0414 


3075 


15.90±0.19 


+0.11 


+1.50±0.84 


16.10±0.16 


+0.31 


+1.49 ±0.74 


HE 0435 


3250 


16.09 ±0.19 


+0.64 


+0.67 ±0.55 


16.37±0.16 


+0.93 


+0.55 ±0.49 


RXJ0911 


2299 


16.23±0.13 


+0.84 


+0.17 ±0.41 


16.39±0.13 


+1.01 


+0.12 ±0.42 


SDSS J0924 


3462 


15.79±0.16 


+0.75 


+0.17 ±0.49 


15.97±0.13 


+0.92 


+0.19 ±0.42 


HE1113 


4438 


15.86±0.18 


+0.74 


+0.05 ± 0.49 


16.03±0.11 


+0.91 


+0.05 ±0.38 


PG1115 


3212 


16.72±0.12 


+1.02 


+0.40 ±0.45 


16.90 ±0.11 


+1.20 


+0.45 ±0.39 


RXJ1131 


5270 


15.55±0.14 


+0.43 


+0.40 ±0.50 


15.80±0.12 


+0.69 


+0.20 ±0.46 


SDSS Jl 138 


2540 


16.01 ±0.19 


+1.44 


+0.41 ± 0.54 


16.26±0.16 


+1.68 


+0.43 ±0.45 


WFIJ2026 


2705 


16.52±0.15 


+1.04 


+0.27 ±0.53 


16.68±0.12 


+1.20 


+0.17 ±0.42 


WFIJ2033 


3285 


16.71 ±0.16 


+1.55 


-0.63 ±0.52 


16.91 ±0.13 


+1.75 


-0.67 ±0.41 



A c is the geometric average of the rest wavelengths of our observations. 



2/3 

a slope that is about half the expected r\/2 oc M B ' H dependence 
(see Equation (J2l)), but with a substantial uncertainty. Perhaps 
more importantly, we find a significant offset from the thin 
disk model in the overall normalization, nearly an order of 
magnitude. The only quasar that we find to be consistent with 
the prediction is MG J0414. 

For compariso n, we ha ve plotted the microlensing size es- 
timates of Morg an et al.| ( |2010) in Figure 10 including their 
wavelength corrections but undoing their inclination correc- 
tions to match our convention and multiplying by a factor of 
2.44 to convert to half-light radii (see Equation Q). They find 
a logarithmic offset between their microlensing sizes and the 
thin disk prediction in the same sense as ours, and about half 
as large. Though we do not have a good explanation for the 
offset between our microlensing size measurements and the 
prediction, the smal ler offset between our sizes and those of 
Morg an et al.| ( |2010| may arise from the difference in analysis 
techniques; analyzing this possibility is the subject of future 
work. 

Increasing the Eddington fraction of the quasars, or de- 
creasing their accretion efficiencies, would help to reconcile 
our microlensing sizes with the predictions of the thin disk 
model, though the adjustments would need to be drastic, be- 



cause of the weak dependence of the size on these values. In- 
creasing the black hole mass estimates systematically by fac- 
tors of ~ 20 or decreasing the mean stellar mass in the lensing 
galaxies by factors of ~ 50 would also bring about agreement, 
but systematic offsets at these levels seem unlikely. 

With one exception, every quasar in our sample displays 
a power-law slope v flatter than the expected 4/3 by 1 .2 to 
2.6cr or more (see Table 15). Although these do not individu- 
ally rule out the thin disk slope at high significance, the 10 to 
1 preponderance of quasars with too-shallow slopes makes a 
convincing case against it. Treating the measurement for each 
quasar as an independent constraint yields an average value 
z/ = 0.17±0.15. Comparing the eleven measured slopes to 
4/3 yields a x 2 /A^of = 6.1. Since a multi-temperature black- 
body accretion disk with a power-law temperature profile has 
?eff oc r~ x l v , applying our results to this model implies that the 
effective temperature is a steeply falling function of radius. It 
is worth noting that nearly all the measured slopes are consis- 
tent with v = 0, i.e., a source with a wavelength-independent 
size, though the existence of such an object seems highly un- 
likely, from both theoretical (e.g., energy conservation) and 
observational (e.g., quasar variability) standpoints. 

MGJ0414, the only quasar that matches the thin disk size 
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o.i 1.0 

Rest Wavelength (microns) 

Figure 9. Medians (squares) and modes (crosses) of the probability distribu- 
tions for the half-light source radius. The error bars correspond to the sum in 
quadratures of the intrinsic width of the distributions and the scatter from the 
Monte Carlo propagation of measurement uncertainties. Black (gray) points 
indicate a logarithmic (linear) prior. The solid line is the prediction for the 
half-light radius (in cm) of the standard thin disk model, for the estimated 
black hole mass. In the case of SDSS J1330, we only have upper limits on 
the size of the disk. 
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Figure 10. Half-light radii at a rest wavelength of 2500 A, plotted against 
black hole mass. Black hole masses estimated using the virial method are 
plotted as red squares; those estimated using bolometric luminosity are plot- 
ted as blue diamonds. The solid line is the prediction of the thin disk model 
(with /ejj = 0.25 and r\ = 0.15), while the dotted line is the best fit to these 
data. The microlensing radii of Morgan et al. 1 20101 are plotted in thin black 
lines for comparison. They have been corrected to 2500 A (assuming that 
r <x A 4 ' 3 ) and converted to half-light radii by multiplying by 2.44 (see Equa- 
tion J5J). The inclination correction has been removed. Their error bars are 
largerbecause they take several sources of systematic error into account. The 
black hole masses have unplotted uncertainties of about 0.4 dex. 



prediction (and the only radio-loud quasar in our sample), is 
also the only quasar with a slope v consistent with 4/3. Al- 
though MG J0414 was the only case where we used empirical 
mid-IR flux ratios in place of model ratios (as we suspected 
the presence of millilensing), it is unlikely that its unique be- 
havior is due to this choice, since mid-IR radios for other 
lensed qua sars (including PG 1 1 15 ) match the models very 
well ( |Chiba et al.|2005||M"inezaki et al.|2009) . 

Two lenses, HE 0230 and WFIJ2033, have microlensing 
size estimates that decrease with wavelength. This result does 
not seem to arise from a gross failure of our analysis method, 
since the flux ratios themselves become more anomalous with 
wavelength, which in most cases implies a smaller source size 
(see Table B). It cannot be due to microlensing variability, 
since the optical/IR measurements are coeval. But it should 
not be concluded that the temperature profile of these quasars 
is in fact inverted, as the errors do not rule out positive slopes. 
It is possible that the sources lie on an unusual region of the 
caustic pattern where the smoothing caused by a larger source 
causes more anomalous flux ratios. These lenses underscore 
the importance of having a large sample of objects when using 
single-epoch measurements. 

8. SUMMARY AND CONCLUSIONS 

We have obtained flux ratios in eight optical and IR filters 
spanning a factor of six in wavelength, for a sample of 12 
strongly lensed quasars. Comparing these ratios to each other 
and to Chandra X-ray flux ratios, both from the literature and 
newly reported here, we see chromatic variations, which we 
attribute to the microlensing of a multicolor accretion disk. 
The standard thin accretion disk model predicts a temperature 
profile that falls as the /? = 3/4 power of radius, implying that 
the half-light radius of the disk is proportional to A 4 / 3 , with an 
overall normalization that depends on the black hole mass and 
(weakly) on the Eddington fraction and accretion efficiency. 
This chromatic dependence combines with the dependence of 
microlensing magnification on source size to produce higher- 
amplitude microlensing variations at blue wavelengths than 
at red wavelengths. For single-epoch observations like ours, 
this usually means that the observed flux ratios at blue wave- 
lengths are more anomalous than at red wavelengths, com- 
pared to the ratios predicted by smooth lens models. 

In addition to formal statistical uncertainties, we have es- 
timated the systematic errors in the flux ratios due to spec- 
tral contamination by broad emission lines, quasar variability 
combined with lensing time delays, and microlensing vari- 
ability combined with delays between measurement epochs. 
We have described a Bayesian method for determining the 
probability distribution of the half-light radius of the source 
in each filter. For quasars without X-ray measurements, we 
are able to place upper limits on the half-light radii, by virtue 
of the departures of the flux ratios from the model predic- 
tions. When X-ray ratios are available, we assume that they 
originate from a very compact region, and we are able to place 
both upper and lower limits on the half-light radius, reasoning 
that the differences between the X-ray and optical ratios arise 
because of the appreciable optical radius. 

The resulting size estimates are larger overall than is pre- 
dicted by the standard thin disk model by nearly an order of 
magnitude, on average. The scaling of the half-light radius 
with mass is consistent with the expected slope, but the scal- 
ing with wavelength is shallower than expected. Although the 
large error bars on the wavelength slopes mean that a single 
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slope measurement would have little weight, the fact that all 
but one of the measured slopes are too shallow commands at- 
tention. Since the scaling of radius with wavelength is the 
reciprocal of the scaling of temperature with radius, our re- 
sults indicate that if we assume a blackbody accretion disk 
with a power-law temperature profile the temperature slope is 
in general steeper than the standard thin disk model predicts. 
However, other explanations for the relative insensitivity of 
radius to wavelength are also possible, though none immedi- 
ately present themselves. 

There is growing evidence from microlensing studies that 
quasar accretion disks are larger at optical wavelengths than 
simple accretion models predict. Ou r result is a co nfirmation 
of earlier, more qualitative w ork by |Pooley et al. ( |2007| l. It 
also corroborates the results of Morgan et al. (2010J), who use 
several lensed quasars to study the mass dependence of the 
accretion disk radius, and find it difficult to reconcile their 
radii with thin disk theory. Our results for the observer-frame 
r'-band size of MGJ0414 and M '-band size of S DSSJ0 924 

20081 and 



microlensing simulation code. J.A.B. and P.L.S. acknowledge 
the support of US NSF grant AST-0607601 . 



are consistent wi th the upper limits of Bate et al. 
|Floyd et"aL] ( |2009| l, respectively. Likewise, our estimates of 
the radius of RXJ1131 at around 4000 A in the rest frame, 
and of HE 0435 at around 6500 A in the obs erved frame, are 
consis tent with those of |Dai et aL] ( |2010| ) and Mosquer a et al.| 
( |2010| ), respectively. 

The picture is less clear regarding the temperature profile 
slope (3. Using time-domain m easurements at several wave- 
lengths, Poindexter et al. (2008) find /3 ~ 0.6 for the quasar 
HE 1 104- 1805, consistent with the thin disk slope of 3/4. 
They point out that a flatter temperature slope (i.e., j3 < 3/4) 
would help to explain the faintness of lensed quasars relative 
to the expected flux from a Shakura-Sunyaev disk of the size 
implied by microlensing. However, our result is just the op- 
posite: the low average value of v implies that the tempera- 
ture slope f3 is greater than 3/4. Other single-lens microlens- 
ing studies have addressed this question as well. |Bate et aTj 
(2008) and Floyd et al. (2009) use single-epoch broadband 
photometry similar to ours to constrain the slope v for the 
quasars MGJ0414 and SDSS J0924, respectively; their lim- 
its are compatible b oth with our r esults and with the thin 
disk model. Mosq uera et al.| ( |2010| l use narrowband imag- 
ing of HE 0435, finding a slope v = 1.3 ±0.3, consistent with 
our slope within the errors. In addition, several studies have 
focused on Q 2237 + 0305, because of the short time scale 
and large amplitude of its microlensing variability, and the 
large amount of monitoring data that are available. Anguita 
et al. ( |2008| ) find a ratio of /- to_g'-band source sizes con- 



sistent with the thin disk model. Eigenbrod et al. ( |2008[ ) use 
spectrophotometric monitoring to constrain v = 1 .2 ± 0.3, and 
Mosquera et al. (2009) use narrowband photometry, finding 
solutions compatible with v = 4/3, as well as with higher val- 
ues. Our result, working with a larger sample size, is the first 
to indicate that v is significantly less than 4/3, so further study 
is clearly needed. We note that the error bars on v are uni- 
formly large for individual lenses. Increasing the sample size, 
as we have done, is a good way to improve this; additional 
strategies include using spectroscopic methods to transform 
line emission from a nuisance into a tool, improving the coor- 
dination between X-ray and optical observing campaigns, and 
moving to multiwavelength, long temporal baseline monitor- 
ing. 



The authors thank Joachim Wambsganss for the use of his 
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Table 5 
Relative Photometry 



Quasar/Filter 



HM a 



HS 



LM 



LS 



HE 0230- 2130 


A 


B 


c 


D 






»' 


OiO.OO 


+0.26 ±0.03 ±0.00 


+0.97 ±0.03 ±0.01 


+2.60 ±0.08 ±0.03 


0.17 




g' 


OiO.OO 


+0.25 ±0.01 ±0.05 


+0.69 ±0.01 ±0.05 


+2.41 ±0.04 ±0.11 


0.14 




r 1 


OiO.OO 


+0.23 ±0.01 ±0.05 


+0.62 ±0.01 ±0.05 


+2.44 ±0.02 ±0.20 


0.17 




i' 


OiO.OO 


+0.17 ±0.01 ±0.00 


+0.58 ±0.01 ±0.01 


+2.37 ±0.05 ±0.36 


0.15 




z' 


OiO.OO 


+0.12±0.01±0.05 


+0.55 ±0.02 ±0.05 


+2.95±0.15±0.57 


0.38 




.1 


OiO.OO 


+0.15 ±0.01 ±0.00 


+0.54 ±0.01 ±0.01 


+2.64 ±0.07 ±0.94 


0.26 




II 


OiO.OO 


+0.12±0.01±0.05 


+0.50 ±0.01 ±0.05 


+2.74±0.12±1.18 


0.30 




K s 


OiO.OO 


+0.07 ±0.01 ±0.05 


+0.50 ±0.01 ±0.05 


+3.48±0.16±1.62 


0.62 


0.5- 


-8keV (ObsID 1642) 


-0.51 ±0.13 ±0.50 


+0.39 ±0.19 ±1.20 


0±0.30 


+0.87±0.18±0.50 


0.60 


MGJ0414 + 0534 


Al 


A2 


B 


C 






T> 


+0.27 ±0.22 ±0.05 


+0.21 ±0.21 ±0.05 


0±0.02 


+0.43 ±0.11 ±0.07 


0.88 




i' 


-0.50 ±0.09 ±0.00 


+0.30±0.16±0.00 


0±0.02 


+1.03 ±0.10 ±0.04 


0.57 




z' 


-0.66 ±0.03 ±0.00 


+0.02 ±0.06 ±0.00 


0±0.02 


+1.02 ±0.05 ±0.04 


0.46 




.1 


-0.93 ±0.01 ±0.00 


-0.41 ±0.01 ±0.00 


0±0.02 


+0.94 ±0.01 ±0.04 


0.31 




II 


-0.96 ±0.01 ±0.05 


-0.72 ±0.01 ±0.05 


0±0.02 


+0.90 ±0.01 ±0.07 


0.21 




K s 


-1.01 ±0.01 ±0.00 


-0.86 ±0.01 ±0.00 


0±0.02 


+0.93 ±0.01 ±0.04 


0.15 


0.5- 


-8keV (ObsID 3919) 


-0.81 ±0.05 ±0.50 


-0.29 ±0.08 ±1.20 


0±0.35 


+0.94 ± 0.05 ± 0.40 


0.36 


HE 0435 -1223 


c 


B 


A 


D 






»' 


+0.93 ±0.02 ±0.01 


+1.01 ±0.02 ±0.01 


0±0.00 


+0.90 ± 0.02 ± 0.02 


0.45 




g' 


+0.73 ±0.01 ±0.05 


+0.77 ±0.01 ±0.05 


0±0.00 


+0.85 ±0.01 ±0.05 


0.35 




r 1 


+0.69 ±0.02 ±0.01 


+0.70 ±0.02 ±0.01 


0±0.00 


+0.83 ±0.03 ±0.02 


0.32 




i' 


+0.56 ±0.01 ±0.05 


+0.62 ±0.01 ±0.05 


0±0.00 


+0.84 ±0.01 ±0.05 


0.28 




z' 


+0.54 ±0.02 ±0.01 


+0.58 ±0.02 ±0.01 


0±0.00 


+0.78 ±0.02 ±0.02 


0.27 




.1 


+0.49 ±0.01 ±0.05 


+0.50 ±0.01 ±0.05 


0±0.00 


+0.69 ±0.01 ±0.05 


0.24 




H 


+0.44 ±0.01 ±0.05 


+0.47 ±0.01 ±0.05 


0±0.00 


+0.67 ±0.01 ±0.05 


0.23 




K s 


+0.36 ±0.01 ±0.01 


+0.38 ±0.01 ±0.01 


0±0.00 


+0.64 ±0.01 ±0.02 


0.19 


0.5- 


-8keV (ObsID 7761) 


+1.06±0.13 


+1.07±0.14±0.35 


0±0.15 


+1.10±0.13±0.20 


0.47 


RXJ091 1+0551 


B 


A 


D 


C 






«' 


+0.37 ± 0.02 ± 0.00 


0±0.00 


+1.64 ±0.01 ±0.09 


+1.31 ±0.02 ±0.00 


0.46 




g' 


+0.26 ±0.01 ±0.10 


0±0.00 


+1.40±0.01±0.13 


+1.18±0.01±0.10 


0.43 




r 1 


+0.28 ±0.01 ±0.05 


0±0.00 


+1.33±0.Ol±0.1O 


+1.14±0.01±0.05 


0.44 




i' 


+0.13 ±0.01 ±0.05 


0±0.00 


+1.20±0.01±0.10 


+0.95 ±0.01 ±0.05 


0.39 




z 1 


+0.24 ±0.01 ±0.00 


0±0.00 


+1.22±0.01±0.09 


+1.00 ±0.01 ±0.00 


0.42 




.1 


+0.06 ± 0.02 ± 0.00 


0±0.00 


+1.09±0.01±0.09 


+0.79 ±0.01 ±0.00 


0.36 




II 


+0.09 ± 0.02 ± 0.00 


0±0.00 


+1.02 ±0.02 ±0.09 


+0.77 ±0.02 ±0.00 


0.39 




K s 


-0.02 ±0.02 ±0.00 


0±0.00 


+0.94 ±0.03 ±0.09 


+0.77 ±0.02 ±0.00 


0.38 


0.5 


-8keV (ObsID 419) 


-0.29 ±0.03 ±0.45 


-1.33±0.17±0.80 


0±0.31 


+1.14±0.37±0.65 


0.98 


SDSSJ0924 + 0219 


A 


D 


B 


C 






«' 


0±0.0() 


+1.66± 0.06 ±0.50 


+1.48 ±0.03 ±0.05 


+2.66 ±0.09 ±0.50 


0.71 




g' 


0±0.00 


+2.76 ±0.03 ±0.05 


+1.53±0.01±0.05 


+2.59 ±0.02 ±0.05 


1.00 




, J 


0±0.00 


+2.96 ±0.03 ±0.00 


+1.43±0.01±0.01 


+2.42 ±0.01 ±0.01 


1.07 




i' 


0±0.00 


+2.96 ±0.03 ±0.05 


+1.36± 0.01 ±0.05 


+2.28 ±0.01 ±0.05 


1.07 




z' 


0±0.00 


+2.85 ±0.04 ±0.00 


+1.36±0.01±0.01 


+2.28 ±0.02 ±0.01 


1.03 




./ 


0±0.00 


+2.23 ±0.06 ±0.05 


+1.16±0.02±0.05 


+1.88±0.04±0.05 


0.81 




II 


0±0.00 


+2.11 ±0.07 ±0.05 


+1.13±0.03±0.05 


+1.65 ±0.04 ±0.05 


0.75 




K s 


± 0.00 


+1.88 ±0.06 ±0.00 


+1.03 ±0.03 ±0.01 


+1.41 ±0.04 ±0.01 


0.68 


0.5- 


-8keV (ObsID 5604) 


-1.26 ±0.23 ±0.30 


+0.87 ±0.54 ±0.70 


0±0.25 


+0.94 ±0.42 ±0.50 


0.78 


HE 1113-0641 


B 


D 


A 


C 






g' 


+0.10 ±0.01 ±0.05 


+0.81 ±0.01 ±0.05 


0±0.00 


+0.43 ±0.01 ±0.05 


0.43 




r 1 


+0.09 ±0.01 ±0.05 


+0.82 ±0.01 ±0.05 


0±0.00 


+0.50 ±0.01 ±0.05 


0.43 




i' 


+0.07 ±0.01 ±0.00 


+0.83 ±0.01 ±0.00 


0±0.00 


+0.45 ±0.01 ±0.00 


0.44 




z' 


+0.06 ±0.01 ±0.00 


+0.86 ±0.01 ±0.00 


0±0.00 


+0.41 ±0.01 ±0.00 


0.46 




.1 


+0.12 ±0.01 ±0.00 


+0.75 ±0.01 ±0.00 


0±0.00 


+0.39 ±0.01 ±0.00 


0.41 




H 


+0.08 ±0.01 ±0.00 


+0.60 ±0.01 ±0.00 


0±0.00 


+0.42 ±0.01 ±0.00 


0.35 




K s 


+0.11 ±0.01 ±0.00 


+0.48 ±0.01 ±0.00 


0±0.00 


+0.55 ±0.01 ±0.00 


0.29 


0.5- 


-8keV (ObsID 7760) 


+0.50 ±0.84 ±0.00 


+0.27 ± 0.47 ± 0.00 


0±0.00 


+1.75±1.19±0.01 


0.52 
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Table 5 
Relative Photometry 



PG 1115+080 


Al 


A2 


c 


B 






«' 


-1.46 ±0.01 ±0.05 


-0.85 ±0.01 ±0.05 


0±0.02 


+0.40 ±0.01 ±0.05 


0.25 




g' 


-1.45 ±0.01 ±0.05 


-1.28±0.01±0.05 


0±0.02 


+0.31 ±0.01 ±0.05 


0.10 




r 1 


-1.41 ±0.01 ±0.00 


-1.23 ±0.01 ±0.00 


0±0.02 


+0.35 ±0.01 ±0.02 


0.11 




i' 


-1.42 ±0.01 ±0.05 


-1.22 ±0.01 ±0.05 


0±0.02 


+0.39 ±0.01 ±0.05 


0.10 




z' 


-1.39±0.01±0.00 


-1.21±0.01±0.00 


0±0.02 


+0.38 ±0.01 ±0.02 


0.10 




.1 


-1.43 ±0.01 ±0.05 


-1.18±0.01±0.05 


0±0.02 


+0.41 ±0.01 ±0.05 


0.11 




II 


-1.43 ±0.02 ±0.05 


-1.17±0.02±0.05 


0±0.02 


+0.34 ±0.02 ±0.05 


0.13 




K s 


-1.40 ±0.03 ±0.00 


-1.16±0.03±0.00 


0±0.02 


+0.48 ±0.05 ±0.02 


0.10 


0.5- 


-8keV (ObsID 7757) 


-2.54 ±0.07 ±0.00 


-2.25 ± 0.09 ± 0.00 


-1.26 ±0.06 ±0.02 


-1.30 ±0.06 ±0.02 


0.34 


RXJ1131-1231 


B 


A 


C 


D 






u' 


-1.29 ±0.01 ±0.00 


-1.74 ± 0.01 ±0.02 


0±0.01 


+0.63 ± 0.02 ± 0.07 


0.74 




g' 


-1.06±0.01±0.05 


-1.52 ±0.01 ±0.05 


0±0.01 


+0.94 ±0.01 ±0.08 


0.62 




r 1 


-1.10 ±0.04 


-1.55 ±0.06 


0±0.01 


+1.00±0.11 


0.54 




i' 


-1.00±0.10 


-1.46±0.12 


0±0.01 


+0.97 ±0.10 


0.58 




z' 


-0.90 ±0.08 


-1.33±0.11 


0±0.01 


+1.00 ±0.08 


0.60 




.1 


-0.78 ±0.17 


-1.21 ±0.23 


0±0.01 


+0.75 ±0.12 


0.74 




II 


-0.64 ±0.1 8 


-1.04 ±0.22 


0±0.01 


+0.65 ±0.14 


0.83 




K s 


-0.51 ±0.12 


-1.12±0.20 


0±0.01 


+0.97 ±0.1 8 


0.77 


0.5- 


-8keV (ObsID 7787) 


-1.26 ±0.08 ±0.00 


-1.89 ±0.08 ±0.02 


0±0.01 


+0.51 ±0.11 ±0.07 


0.78 


0.5- 


-8keV (ObsID 7789) 


-1.16 ±0.08 ±0.00 


-1.87 ±0.07 ±0.02 


0±0.01 


+1.07 ±0.12 ±0.07 


0.63 


SDSSJ1138 + 0314 


A 


D 


c 


B 






u' 


0±0.01 


+1.98 ±0.05 ±0.05 


+1.70±0.03±0.05 


+1.35 ±0.03 ±0.05 


0.72 




g' 


0±0.01 


+1.16±0.01±0.10 


+1.28±0.01±0.10 


+1.33±0.01±0.10 


0.42 




r 1 


0±0.01 


+1.28±0.01±0.05 


+1.27±0.01±0.05 


+1.38±0.01±0.05 


0.45 




i' 


0±0.01 


+1.36±0.01±0.01 


+1.28 ±0.01 ±0.00 


+1.37±0.01±0.01 


0.47 




z' 


0±0.01 


+1.17±0.01±0.05 


+1.19±0.01±0.05 


+1.21 ±0.01 ±0.05 


0.41 




.1 


0±0.01 


+1.33 ±0.02± 0.01 


+1.16±0.01±0.00 


+1.39 ±0.02 ±0.01 


0.45 




II 


0±0.01 


+1.26 ±0.03 ±0.05 


+1.06±0.02±0.05 


+1.33 ±0.03 ±0.05 


0.42 




K s 


0±0.01 


+0.83 ±0.03 ±0.05 


+0.89 ±0.03 ±0.05 


+1.23 ±0.04 ±0.05 


0.27 


0.5- 


-8keV (ObsID 7759) 


-1.26±0.35±0.01 


-0.29 ±0.47 ±0.01 


0±0.00 


+0.00 ±0.46 ±0.01 


0.37 


WFIJ2026-4536 


Al 


A2 


B 


C 






«' 


-0.61 ±0.01 ±0.00 


-1.30 ± 0.01 ±0.00 


0±0.02 


+0.18 ±0.01 ±0.01 


0.34 




g' 


-0.99 ±0.01 ±0.05 


-1.27 ±0.01 ±0.05 


0±0.02 


+0.22 ±0.01 ±0.05 


0.20 




r 1 


-1.25 ±0.01 ±0.05 


-1.21 ±0.01 ±0.05 


0±0.02 


+0.22 ±0.01 ±0.05 


0.12 




i> 


-1.35±0.01±0.00 


-1.17±0.01±0.00 


0±0.02 


+0.23 ±0.01 ±0.01 


0.12 




z' 


-1.38±0.01±0.05 


-1.13 ±0.01 ±0.05 


0±0.02 


+0.25 ±0.01 ±0.05 


0.12 




.1 


-1.53±0.01±0.00 


-1.11±0.01±0.00 


0±0.02 


+0.29 ±0.01 ±0.01 


0.18 




II 


-1.49 ±0.01 ±0.07 


-0.90 ±0.01 ±0.07 


0±0.02 


+0.30 ±0.02 ±0.05 


0.19 




K s 


-1.45 ±0.01 ±0.05 


-0.83 ±0.02 ±0.05 


0±0.02 


+0.26 ±0.02 ±0.05 


0.19 


0.5- 


-8keV (ObsID 7758) 


-1.91±0.34±0.20 


-0.75±0.71±1.55 


0±0.10 


+1.00±0.33±0.15 


0.58 


WFIJ2033-4723 


Al 


A2 


B 


C 






u' 


-0.40 ±0.02 ±0.00 


+0.07 ± 0.03 ± 0.00 


0±0.03 


+0.70 ± 0.03 ± 0.03 


0.07 




g' 


-0.53 ±0.01 ±0.05 


+0.14±0.01±0.05 


0±0.03 


+0.48 ±0.01 ±0.06 


0.06 




r 1 


-0.52 ±0.01 ±0.00 


+0.05 ±0.01 ±0.00 


0±0.03 


+0.38 ±0.01 ±0.03 


0.06 




i' 


-0.55 ±0.01 ±0.05 


-0.05 ±0.01 ±0.05 


0±0.03 


+0.23 ±0.01 ±0.06 


0.10 




z' 


-0.53 ±0.02 ±0.00 


-0.05 ±0.02 ±0.00 


0±0.03 


+0.28 ±0.02 ±0.03 


0.08 




.1 


-0.60 ±0.01 ±0.05 


-0.12±0.01±0.05 


0±0.03 


+0.15 ±0.01 ±0.06 


0.12 




II 


-0.58 ±0.01 ±0.05 


-0.14±0.01±0.05 


0±0.03 


+0.13 ±0.01 ±0.06 


0.13 




K s 


-0.58 ±0.01 ±0.00 


-0.10±0.02±0.00 


0±0.03 


+0.10 ±0.02 ±0.03 


0.15 


0.5 


-8keV (ObsID 5603) 


+0.15±0.18±0.25 


+0.00 ±0.22 ±0.35 


0±0.20 


+0.49 ±0.17 ±0.25 


0.28 



Note. — All values are in magnitudes. Formal statistical uncertainties are reported first, followed by estimated systematic uncertainties. 
Where there are no formal uncertainties, only the systematics are reported. 
a HS: highly magnified saddle point; HM: highly magnified minimum; LS: less magnified saddle point; LM: less magnified minimum. 

Standard deviation of the difference of the four image magnitudes and the lens model prediction in magnitudes. This quantity serves as a 
quick estimate of the degree of anomaly in the flux ratios. 
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Table 7 












Half-light Radii 








Quasar /Wavelength 


logOi/ 2 /cm), Loj 


; Prior 


log(>i/2/ 


cm), Linear Prior 


(fim) 


Median 


Mode 


Error 


Median 


Mode 


Error 


HE 0230-2130 


0.11 


16.66 


16.78 


0.26 


16.75 


16.83 


0.22 


0.15 


16.80 


16.89 


0.34 


16.92 


16.96 


0.29 


0.20 


16.70 


16.81 


0.34 


16.82 


16.87 


0.32 


0.24 


16.65 


16.75 


0.34 


16.77 


16.83 


0.33 


0.29 


16.43 


16.57 


0.34 


16.53 


16.61 


0.30 


0.40 


16.44 


16.57 


0.36 


16.55 


16.62 


0.35 


0.52 


16.35 


16.35 


0.38 


16.47 


16.60 


0.36 


0.68 


16.29 


16.32 


0.42 


16.40 


16.59 


0.40 


MGJ0414 + 0534 


0.17 


15.38 


15.50 


0.34 


15.52 


15.62 


0.31 


0.21 


15.68 


15.81 


0.45 


15.85 


15.95 


0.36 


0.25 


15.77 


15.92 


0.41 


15.94 


16.05 


0.38 


0.34 


15.96 


16.12 


0.40 


16.10 


16.21 


0.34 


0.45 


16.13 


16.29 


0.43 


16.30 


16.40 


0.35 


0.59 


16.23 


16.36 


0.38 


16.37 


16.46 


0.36 


HE 0435 -1223 


0.13 


15.71 


15.82 


0.42 


15.97 


16.03 


0.38 


0.18 


15.89 


15.98 


0.42 


16.15 


16.24 


0.37 


0.23 


15.93 


16.04 


0.38 


16.17 


16.27 


0.32 


0.29 


16.03 


16.16 


0.37 


16.27 


16.37 


0.32 


0.34 


16.09 


16.26 


0.36 


16.31 


16.39 


0.30 


0.46 


16.17 


16.35 


0.36 


16.39 


16.47 


0.31 


0.61 


16.20 


16.38 


0.36 


16.41 


16.49 


0.31 


0.80 


16.28 


16.43 


0.34 


16.47 


16.52 


0.33 


RXJ091 1+0551 


0.09 


16.14 


16.20 


0.32 


16.30 


16.29 


0.33 


0.13 


16.10 


16.17 


0.31 


16.25 


16.26 


0.31 


0.17 


16.13 


16.19 


0.32 


16.29 


16.27 


0.34 


0.20 


16.16 


16.21 


0.30 


16.32 


16.31 


0.34 


0.24 


16.19 


16.24 


0.30 


16.35 


16.34 


0.33 


0.33 


16.26 


16.32 


0.29 


16.40 


16.40 


0.32 


0.43 


16.17 


16.25 


0.30 


16.31 


16.33 


0.32 


0.57 


16.25 


16.32 


0.27 


16.37 


16.38 


0.28 


SDSSJ0924 + 0219 


0.14 


15.78 


15.92 


0.35 


15.93 


16.05 


0.31 


0.19 


15.65 


15.80 


0.33 


15.78 


15.91 


0.28 


0.25 


15.66 


15.83 


0.32 


15.78 


15.92 


0.26 


0.31 


15.68 


15.85 


0.32 


15.79 


15.93 


0.26 


0.36 


15.68 


15.84 


0.32 


15.80 


15.93 


0.26 


0.50 


15.75 


15.87 


0.37 


15.90 


15.98 


0.31 


0.65 


15.79 


15.91 


0.39 


15.94 


16.02 


0.32 


0.86 


15.84 


15.97 


0.34 


16.00 


16.08 


0.30 


HE 11 13 -0641 


0.21 


15.78 


15.96 


0.31 


15.93 


16.04 


0.23 


0.28 


15.84 


16.02 


0.33 


15.99 


16.11 


0.23 


0.34 


15.84 


16.02 


0.33 


15.99 


16.10 


0.23 


0.41 


15.78 


15.94 


0.30 


15.91 


16.03 


0.23 


0.56 


15.75 


15.93 


0.31 


15.90 


16.01 


0.23 


0.74 


15.83 


16.01 


0.31 


15.97 


16.09 


0.22 


0.97 


15.87 


16.05 


0.34 


16.02 


16.14 


0.23 


PG 1115 + 080 


0.13 


16.34 


16.49 


0.32 


16.48 


16.57 


0.25 


0.18 


16.65 


16.72 


0.37 


16.78 


16.81 


0.35 


0.23 


16.69 


16.79 


0.30 


16.81 


16.85 


0.28 


0.28 


16.76 


16.82 


0.35 


16.89 


16.89 


0.32 


0.33 


16.74 


16.82 


0.30 


16.86 


16.89 


0.29 


0.46 


16.74 


16.81 


0.34 


16.86 


16.86 


0.31 


0.61 


16.63 


16.73 


0.32 


16.75 


16.79 


0.27 


0.79 


16.84 


16.88 


0.32 


16.96 


17.11 


0.29 
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Table 7 


















Half-light Radii 








RXJ1131-1231 


0.22 


15.35 


15.45 


0.35 


15.62 


15.73 


0.40 


0.29 


15.41 


15.54 


0.47 


15.69 


15.81 


0.35 


0.38 


15.36 


15.49 


0.47 


15.66 


15.78 


0.35 


0.46 


15.48 


15.62 


0.49 


15.75 


15.86 


0.35 


0.55 


15.57 


15.70 


0.42 


15.81 


15.90 


0.31 


0.75 


15.59 


15.69 


0.32 


15.77 


15.87 


0.26 


1.00 


15.58 


15.67 


0.29 


15.75 


15.84 


0.25 


1.30 


15.66 


15.77 


0.32 


15.85 


15.94 


0.28 


SDSSJ1138 + 0314 


0.10 


15.78 


15.88 


0.33 


15.95 


16.06 


0.30 


0.14 


15.92 


16.04 


0.43 


16.13 


16.19 


0.34 


0.18 


15.92 


16.04 


0.39 


16.12 


16.19 


0.31 


0.22 


15.92 


16.04 


0.37 


16.12 


16.20 


0.31 


0.26 


15.96 


16.09 


0.51 


16.17 


16.24 


0.40 


0.36 


15.96 


16.10 


0.38 


16.17 


16.25 


0.31 


0.48 


16.01 


16.16 


0.40 


16.22 


16.31 


0.33 


0.63 


16.23 


16.35 


0.45 


16.43 


16.47 


0.37 


SDSSJ1330+1810 a 


0.15 


15.96 












16.22 










0.20 


15.98 












16.25 










0.26 


16.15 












16.42 










0.32 


16.11 












16.38 










0.38 


16.05 












16.33 










0.52 


16.24 












16.49 










0.69 


16.23 












16.49 










0.90 


16.23 












16.55 










WFIJ2026-4536 


0.11 


15.73 


15.83 


0.52 


16.07 


16.08 


0.44 


0.15 


16.20 


16.31 


0.44 


16.46 


16.56 


0.37 


0.20 


16.47 


16.56 


0.32 


16.61 


16.67 


0.27 


0.24 


16.79 


16.88 


0.34 


16.91 


16.98 


0.29 


0.28 


16.52 


16.63 


0.30 


16.65 


16.73 


0.25 


0.39 


16.36 


16.53 


0.32 


16.51 


16.61 


0.22 


0.51 


16.34 


16.52 


0.47 


16.56 


16.69 


0.29 


0.67 


16.27 


16.44 


0.47 


16.53 


16.67 


0.30 


WFIJ2033-4723 


0.14 


16.86 


16.98 


0.36 


16.97 


17.04 


0.28 


0.18 


16.89 


16.96 


0.40 


17.02 


17.04 


0.35 


0.24 


16.84 


16.91 


0.37 


16.96 


17.01 


0.31 


0.29 


16.58 


16.63 


0.38 


16.71 


16.81 


0.31 


0.34 


16.65 


16.69 


0.38 


16.78 


16.87 


0.30 


0.47 


16.51 


16.58 


0.39 


16.66 


16.75 


0.30 


0.62 


16.47 


16.54 


0.37 


16.61 


16.70 


0.29 


0.81 


16.44 


16.52 


0.36 


16.58 


16.69 


0.28 



1 90% upper limits on the half-light radius for SDSS J 1 330. Because our priors cut off below a certain source size, these limits are conservative. 



